Frustration and entanglement in the t2g spin— orbital model on a triangular lattice: 

valence— bond and generalized liquid states 



Bruce Normand 

Departement de Physique, Universite de Fribourg, CH-1700 Fribourg, Switzerland 
Theoretische Physik, ETH-Honggerberg, CH-8093 Zurich, Switzerland 

Andrzej M. Oles 

Marian Smoluchowski Institute of Physics, Jagellonian University, Reymonta 4, PL-30059 Krakow, Poland 
Max- Planck- Institut fiir Festkorperforschung, Heisenbergstrasse 1, D-70569 Stuttgart, Germany 

(Dated: May 7, 2008) 

We consider the spin-orbital model for a magnetic system with singly occupied but triply degen- 
erate t2g orbitals coupled into a planar, triangular lattice, as would be exemplified by NaTi02. We 
investigate the ground states of the model for interactions which interpolate between the limits of 
pure superexchange and purely direct exchange interactions. By considering ordered and dimerized 
states at the mean-field level, and by interpreting the results from exact diagonalization calculations 
on selected finite systems, we demonstrate that orbital interactions are always frustrated, and that 
orbital correlations are dictated by the spin state, manifesting an intrinsic entanglement of these 
degrees of freedom. In the absence of Hund coupling, the ground state changes from a highly res- 
onating, dimer-based, symmetry-restored spin and orbital liquid phase, to one based on completely 
static, spin-singlet valence bonds. The generic properties of frustration and entanglement survive 
even when spins and orbitals are nominally decoupled in the ferromagnetic phases stabilized by a 
strong Hund coupling. By considering the same model on other lattices, we discuss the extent to 
which frustration is attributable separately to geometry and to interaction effects. 



PACS numbers: 71.10.Fd, 74.25. Ha, 74.72.-h, 75.30.Et 
I. INTRODUCTION 

Frustration in magnetic systems may be of geometri- 
cal origin, or may arise due to competing exchange in- 
teractions, or indeed both.l For quantum spins, frustra- 
tion acts to enhance the effects of quantum fluctuations, 
leading to a number of different types of magnetically 
disordered state, among which some of the more familiar 
are static and resonating valence-bond (VB) phases. A 
further form of solution in systems with frustrated spin 
interactions is the emergence of novel ordered states from 
a highly degenerate manifold of disordered states, and the 
mechanism for their stabilization has become known sim- 
ply as "order-by-disorder" Ja^ Many materials are now 
known whose physical properties could be understood 
only by employing microscopic models with frustrated 
spin interactions in which some of these theoretical con- 
cepts operate. 

A different and still richer situation occurs in the class 
of transition-metal oxides or fluorides with partly filled 
2id orbitals and near-degeneracy of active orbital degrees 
of freedom. In undoped systems, large Coulomb interac- 
tions on the transition-metal ions localize the electrons, 
and the low-energy physics is that of a Mott (or charge- 
transfer^) insulator. Their magnetic properties are de- 
scribed by superexchange spin-orbital models, derived 
directly from the real electronic structure and contain- 
ing linearly independent but strongly coupled spin and 
orbital operators,^ Such models emerge from the charge 
excitations which involve various multiplet states,— in 
which ferromagnetic (FM) and antiferromagnetic (AF) 



interactions, as well the tendencies towards ferro-orbital 
(FO) and alternating orbital (AO) order, compete with 
each other. This leads to a profound, intrinsic frustration 
of spin-orbital exchange interactions, which occurs even 
in case of only nearest-neighbor interactions for lattices 
with unfrustrated geometry, such as the square and cu- 
bic lattices.^ The underlying physics is formulated in the 
Goodenough-Kanamori rules fS. which imply that the two 
types of order are complementary in typical situations: 
AO order favors a FM state while FO order coexists with 
AF spin order. Only recently have exceptions to these 
rules been noticedj^ and the search for such exceptions, 
and thus for more complex types of spin-orbital order or 
disorder, have become the topic of much active research. 

A case study for frustration in coupled spin-orbital 
systems is provided by the one-dimensional (ID) SU(4) 
modeliiS One expects a priori no frustration in one di- 
mension and with only nearest -neighbor interactions. 
However, spin and orbital interactions, the latter for- 
mulated in terms of pseudospin operators, appear on a 
completely symmetrical footing for every bond, and favor 
respectively AF and AO ordering tendencies, which com- 
pete with each other. In fact a low-energy but magnet- 
ically disordered spin state also frustrates the analogous 
pseudospin-disordered state, and conversely. This com- 
petition results in strong, combined spin-orbital quan- 
tum fluctuations which make it impossible to separate 
the two subsystems, and it is necessary to treat explicitly 
entangled spin-pseudospin states.—!^ While in one sense 
this may be considered as a textbook example of frustra- 
tion and entanglement, the symmetry of the entangled 
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sectors is so high that joint spin-pseudospin operators are 
as fundamental as the separate spin and pseudospin oper- 
ators, forming parts of a larger group of elementary (and 
disentangled) generators. The fact that the ID SU(4) 
model is exactly solvable also results in fundamental sym- 
metries between the intcrsite correlation functions for the 
spin and orbital (and spin-orbital) sectors^ We return 
below to a more detailed discussion of entanglement and 
its consequences. Although indicative of the rich under- 
lying physics (indeed, unconventional behavior has been 
identified for the SU(4) Hamiltonian on the triangular 
lattice^Siii^) the implications of this model are rather lim- 
ited because it does not correspond to the structure of 
superexchange interactions in real correlated materials. 

Realistic superexchange models for perovskite 
transition-metal oxides with orbital degrees of freedom 
have been known for more than three decades,^ but 
the intrinsic frustrating effects of spin-orbital interac- 
tions have been investigated only in recent yearsi^ii^ 
A primary reason for this delay was the complexity 
of the models and the related quantum phenomena, 
which require advanced theoretical methods beyond a 
straightforward mean-field theory. The structure of 
spin-orbital superexchange involves interactions between 
SU(2)-symmetric spins {Si, Sj} on two nearest-neighbor 
transition-metal ions each coupled to orbital op- 

erators {Ti,Tj} which obey only much lower symmetry 
(at most cubic for a cubic lattice), and its general form 
isi 

^^ = ^E {j'H^^ ■ + ^^S'} ■ (1-1) 

The energy scale J is determined (Sec. II) by the in- 
teraction terms and effective hopping matrix elements 
between pairs of directional Cg orbitals [{dda) element] 

or t2g orbitals [{ddn) element] The orbital operators J^j ^ 

and K^j'' specify the orbitals on each bond (ij) || 7, which 

participate in ^ d"^^dj~^ virtual excitations, and 
thus have the symmetry of the lattice. The form of the 
orbital operators depends on the valence n, on the type 
(eg or t2g) of the orbitals and, crucially, on the bond di- 
rection in real space. It is clear from Eq. (|1.1|) that indi- 
vidual terms in the Hamiltonian Tij can be minimized for 
particularly chosen spin and orbital configurations,^ but 
in general the structure of the orbital operators ensures 
a competition between the diff'erent bonds. 

This directional nature is the microscopic origin of the 
intrinsic frustration mentioned above, which is present 
even in the absence of geometrical frustration. Both 
spin and orbital interactions are frustrated, making long- 
range order more difficult to realize in cither sector, and 
enhancing the effects of quantum fiuctuations. Quite 
generally, because insufficient (potential) energy is avail- 
able from spin or orbital order, instead the system is 
driven to gain (kinetic) energy from resonance processes, 
promoting phases with short-range dynamical correla- 
tions and leading naturally to spin and/or orbital dis- 



order. Disordering tendencies are particularly strong in 
highly symmetric systems, which for crystalline materials 
means cubic and hexagonal structures. Among possible 
magnetically disordered phases for spin systems, tenden- 
cies towards dimer formation are common in the regime 
of predominantly AF spin interactions, and new phases 
with VB correlations occur. This type of physics was 
discussed first for eg orbitals on the cubic lattice,^ and, 
in the context of BaVSa, for one version of the problem 
of t2g orbitals on a triangular latticeii^ The same generic 
behavior has since been found for t2g orbitals on the cubic 
latticcjii Cg-orbital systems on the triangular latticer^ii^ 
and for t2g orbitals in the pyrochlore geometryi22i2i By 
analogy with spin liquids, the orbital-liquid phased has 
been introduced for systems with both e^^'^^ and t2g^^^ 
orbital degrees of freedom. The orbital liquid is a phase 
in which strong orbital fluctuations restore the symmetry 
of the orbital sector, in the sense that the instantaneous 
orbital state of any site is pure, but the time average 
is a uniform occupation of all available orbital states. 
We note that in the discussion of orbital liquids in t2g 
systemsr^i^ it was argued that the spin sector would be 
ordered. To date little is known concerning the behavior 
of orbital correlations in an orbital liquid, the possible in- 
stabilities of the orbital liquid towards dimerized or VB 
phases, or its interplay with lattice degrees of freedom. 

One possible mechanism for the formation of an orbital 
liquid state is the positional resonance of VBs. There has 
been considerable recent discussion of spin-orbital mod- 
els in the continuing search for a realistic system real- 
izing such a resonating VB (RVB) stater^ including in 
a number of the references cited in the previous para- 
graph. While the RVB state was first proposed for the 
S = \ Heisenberg model on a triangular lattice^! exten- 
sive analysis of spin-only models has not yet revealed 
a convincing candidate system, although the nearest- 
neighbor dimer basis has been shown to deliver a very 
good description of the low-energy sector for the S — ^ 
Heisenberg model on a kagome lattice.^^ To date, the 
only rigorous proof for RVB states has been obtained in 
rather idealized quantum dimer models (QDMs)^ most 
notably on the triangular latticci^ The insight gained 
from this type of study can, however, be usecU^ to for- 
mulate some qualitative criteria for the emergence of an 
RVB ground state. These combine energetic and topo- 
logical requirements, both of which are essential: the 
energetics of the system must establish a proclivity for 
dimer formation, a high quasi-degeneracy of basis states 
in the candidate ground manifold, and additional energy 
gains from dimer resonance; exact degeneracy between 
topological sectors (determined by a non-local order pa- 
rameter related to winding of wave functions around the 
system) is a prerequisite to remove the competing possi- 
bility of a "solid" phase with dimer, plaquette or other 
"crystalline" order. ^® 

We comment here that the "problem" of frustration, 
and the resulting highly degenerate manifolds of states 
which may promote resonance phenomena, is often solved 
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FIG. 1: (Color online) Structure of the transition-metal oxide 
with edge-sharing octahedra realized for NaTi02: (a) frag- 
ment of crystal structure, with Ti and Na ions shown respec- 
tively by black and green (grey) circles separated by O ions 
(open circles); (b) titanium (111) plane with adjacent oxy- 
gen layers, showing each Ti^^ ion coordinated by six oxygen 
atoms (open circles). The directions of the Ti-Ti bonds are 
labeled as XY, YZ, and ZX, corresponding to the plane 
spanned by the connecting Ti-O bonds. This figure is repro- 
duced from Ref. [3^. where it served to explain the structure 
of LiNi02. 



by interactions with the lattice. Lattice deformations act 
to lift degeneracies and to stabilize particular patterns 
of spin and orbital order, the most familiar situation 
being that in colossal-magnetoresistance manganites;2£ 
The same physics is also dominant in a number of spinels, 
where electron-lattice interactions are responsible both 
for the Verwey transition in magnetite'^'^ and for t2g or- 
bital order below it, as well as for inducing the Peierls 
state in CuIr2S4 and MgTi204."^^ Similar phenomena are 
also expected^i to play a role in NaTi02. Here, however, 
we will not introduce a coupling to phonon degrees of 
freedom, and focus only on purely electronic interactions 
whose frustration is not quenched by the lattice. 

The spin-orbital interactions on a triangular lattice 
are particularly intriguing. This lattice occurs for edge- 
sharing MOg octahedra in structures such as NaNi02 or 
LiNi02, where the consecutive (111) planes of Ni'^^ ions 
are well separated. These two e^-electron systems be- 
have quite differently: while NaNi02 undergoes a cooper- 
ative Jahn-Teller structural transition followed by a mag- 
netic transition at low temperatures (Tn = 20 K), both 
transitions are absent in LiNi02i^ Possible reasons for 
this remarkable contrast were discussed in Ref. ^33, where 
the authors noted in particular that realistic spin-orbital 
superexchange neither has an SU(2)(g)SU(2) structure, 
nor can it ever be reduced only to the consideration of 
FM spin terms These studies showed in addition that 
LiNi02 is not a spin-orbital liquid, and that the rea- 
sons for the observed disordered state are subtle, as spins 
and orbitals are thought likely to order in a strictly two- 
dimensional (2D) spin-orbital modeli^ 

The possibilities offered for exotic phases in this type of 
model and geometry motivate the investigation of a real- 
istic spin-orbital model with active t2g orbitals, focusing 
first on 3d^ electronic configurations. The threefold de- 



generacy of the orbitals is maintained, although, as noted 
above, this condition may be hard to maintain in real 
materials at low temperatures. A material which should 
exemplify this system is NaTi02 (Fig. [T|) , which is com- 
posed of Ti'^^ ions in configuration, but has to date 
had rather limited experimental^i^S and theoretical^ 
attention. Considerably more familiar is the set of tri- 
angular cobaltates best known for superconductivity in 
Naa;Co02: here the 00"*+ ions have configuration and 
are expected to be analogous to the case by particle- 
hole symmetry. The effects of doping have recently been 
removed by the synthesis of the insulating end-member 
Co02i^ Another system for which the same spin-orbital 
model could be applied is Sr2V04, where the V''+ ions 
occupy the sites of a square lattice<^ 

The model with hopping processes of pure superex- 
change type was considered in the context of doped cobal- 
tates by Koshibae and Maekawa^S These authors noted 
that, like the cubic system, two t2g orbitals are active 
for each bond direction in the triangular lattice, but that 
the superexchange interactions are very different from the 
cubic case because the effective hopping interchanges the 
active orbitals. Here we focus only on insulating systems, 
whose entire low-energy physics is described by a spin- 
orbital model. In addition to superexchange processes 
mediated by the oxygen ions, on the triangular lattice 
it is possible to have direct-exchange interactions, which 
result from charge excitations due to direct d — d hop- 
ping between those t2g orbitals which do not participate 
in the superexchange. The ratio of these two types of 
interaction (a, defined in Sec. II) is a key parameter of 
the model. Further, in transition-metal ions'* the coef- 
ficients of the different microscopic processes depend on 
the Hund exchange Jh arising from the multiplet struc- 
ture of the excited intermediate d^ state, '^^ and we intro- 
duce 



U '' 



(1.2) 



as the second parameter of the model. The aim of this 
investigation is to establish the general properties of the 
phase diagram in the (a, r\) plane. 

We conclude our introductory remarks by returning to 
the question of entanglement. In the analysis to follow 
we will show that the presence of conflicting ordering ten- 
dencies driven by different components of the frustrated 
intersite interactions can be related to the entanglement 
of spin and orbital interactions. By "entanglement" we 
mean that the correlations in the ground state involve 
simultaneous fluctuations of the spin and orbital com- 
ponents of the wave function which cannot be factor- 
ized. We will introduce an intersite spin-orbital corre- 
lation function to identify and quantify this type of en- 
tanglement in different regimes of the phase diagram. It 
has been shownS. that such spin-orbital entanglement is 
present in cubic titanates or vanadates for small values 
of the Hund exchange r\. Here we will find entanglement 
to be a generic feature of the model for all exchange in- 



4 



teractions, even in the absence of dimer resonance, and 
that only the FM regime at sufficiently high rj, which 
is fully factorizable, provides a counterpoint where the 
entanglement vanishes. 

The paper is organized as follows. In Sec.|TT]we derive 
the spin-orbital model for magnetic ions with the elec- 
tronic configuration (Ti'^+ or V^+) on a triangular lattice. 
The derivation proceeds from the degenerate Hubbard 
model, and the resulting Hamiltonian contains both su- 
perexchange and direct exchange interactions. We begin 
our analysis of the model, which covers the full range 
of physical parameters, in Sec. IIIII by considering pat- 
terns of long-ranged spin and orbital order representa- 
tive of all competitive possibilities. These states compete 
with magnetically or orbitally disordered phases domi- 
nated by VB correlations on the bonds, which are in- 
vestigated in Sec. IIVI The analysis suggests strongly 
that all long-range order is indeed destabilized by quan- 
tum fluctuations, leading over much of the phase dia- 
gram to liquid phases based on fluctuating dimers, with 
spin correlations of only the shortest range. In Sec. |V] 
we present the results of exact diagonalization calcula- 
tions performed for small clusters with three, four, and 
six bonds, which reinforce these conclusions and provide 
detailed information about the local physical processes 
leading to the dominance of resonating dimer phases. In 
each of Sees. IIIII ITVl and |Vl we conclude with a short 
summary of the primary results, and the reader who is 
more interested in an overview, rather than in detailed 
energetic comparisons and actual correlation functions 
for the different phases, may wish to read only these. 
Some insight into the competition and collaboration be- 
tween frustration effects of different origin can be ob- 
tained by varying the geometry of the system, and Sec.lVll 
discusses the properties of the model on related lattices. 
A discussion and concluding summary are presented in 
Sec.lvni 



II. SPIN-ORBITAL MODEL 

A. Hubbard model for t2g electrons 

We consider the spin-orbital model on the triangular 
lattice which follows from the degenerate Hubbard-like 
model for t2g electrons. It contains the electron kinetic 
energy and electronic interactions for transition-metal 
ions arranged on the (111) planes of a compound with 
local cubic symmetry and with the ionic configuration, 
and as such is applicable to Ti'^+ or [Fig. 1(a)]. The 
kinetic energy is given by 

^* = - E + (2-1) 

(y>ll7,A"',o- 

where dl^^ are creation operators for an electron with 
spin a =t, i and orbital "color" fi at site i, and the sum is 




(b) (c) 

FIG. 2: (Color online) (a) Schematic representation of the 
hopping processes in Eq. (|2.ip which contribute to magnetic 
interactions on a representative bond {ij) along the c-axis 
in the triangular lattice. The t2g orbitals are represented by 
different colors (greyscale intensities). Superexchange pro- 
cesses involve O 2pz orbitals (violet), and couple pairs of a 
and b orbitals (red, green) with effective hopping elements t, 
interchanging their orbital color. Direct exchange couples c 
orbitals (blue) with hopping strength t' . (b) Pairs of t2g or- 
bitals active in superexchange and (c) single orbitals active in 
direct exchange; horizontal bonds correspond to the situation 
depicted in panel (a). 



made over all the bonds {ij)\\"f spanning the three direc- 
tions, J = a,b, c, of the triangular lattice. This notation 
is adopted from the situation encountered in a cubic array 
of magnetic ions, where only two of the three orbitals 

are active on any one bond (ij), and contribute 1^/2) to the 
kinetic energy, while the third lies in the plane perpendic- 
ular to the 7 axis and thus hopping processes involving 
the 2p7r oxygen orbitals is forbidden by symmetryj ^^i^^ 
We introduce the labels a = yz, b = xz, and c= xy also 
for the three orbital colors, and in the figures to follow 
their respective spectral colors will be red, green, and 
blue. 

For the triangular lattice formed by the ions on the 
(111) planes of transition-metal oxides (Fig. [J) it is also 
the case that only two t2g orbitals participate in (superex- 
change) hopping processes via the oxygen sites. However, 
unlike the cubic lattice, where the orbital color is con- 
served, here any one active orbital color is exchanged for 
the other one [Fig. [H^a)]. Using the same convention, 
that each direction in the triangular lattice is labeled by 
its inactive orbital color 7 = a,6, c, the hopping ele- 
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ments for a bond oriented (for example) along the c-axis 



in Eq. 1^ are t^^^ = t^^J = t, while tia^ = = 0. In 



addition, and also in contrast to the cubic system, for the 
triangular geometry a direct hopping from one c orbital 
to the other, i.e. without involving the oxygen orbitals, is 
also permitted on this bond (Fig. [2]), and this element is 
denoted by t' = t^cc ■ We will also refer to these hopping 
processes as off-diagonal and diagonal. We stress that 
while the lattice structure of magnetic ions is triangular, 
the system under consideration retains local cubic sym- 
metry of the metal-oxygen octahedra, which is crucial 
to ensure that the degeneracy of the three t2g orbitals is 
preserved. 

The electron-electron interactions are described by the 
on-site terms^ 



which gives Js = Jcos^O, Jm = J sin 9 cos 9, and Jd = 
Jsin^ 9; J is the energy unit, which specifies respectively 
the superexchange (J = Jg) and direct-exchange (J — 
Jd) constants in the two limits a ~ and a = 1. The 
Hamiltonian 

H = J {(1 - a) + v/(l - ah 'Hm + a Hd] (2.7) 

consists of three terms which follow from the processes 
described by the exchange elements in Eqs. (|2.4p . each of 
which contains contributions from both high- and low- 
spin excitations. 



B. Superexchange 



where U and Jh represent respectively the intraorbital 
Coulomb and on-site Hund exchange interactions. Each 
pair of orbitals {/z, v} is included only once in the inter- 
action terms. The Hamiltonian (|2.2p describes rigorously 
the multiplet structure of (P ions within the subspace, 
and is rotationally invariant in the orbital space4^ 

When the Coulomb interaction is large compared with 
the hopping elements {U ^ t,t'), the system is a Mott 
insulator with one d electron per site in the t2g orbitals, 
whence the local constraint in the strongly correlated 
regime is 



nib 



1, 



(2.3) 



where 
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di- 



stricted space rii^ = 0, 1. 



The operators act in the re- 
The low-energy Hamiltonian 



may be obtained by second-order perturbation theory, 
and consists of a superposition of terms which follow from 
virtual d\d}j ^ dfd^ excitations. Because each hopping 
process may be of either off-diagonal (t) [Fig. [Hb)] or 
diagonal {t') type [Fig. HJc)], the Hamiltonian consists 
of several contributions which are proportional to three 
coupling constants. 
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(2.4) 



These represent in turn the superexchange term, the di- 
rect exchange term, and mixed interactions which arise 
from one diagonal and one off-diagonal hopping process. 

We choose to parameterize the Hamiltonian by the sin- 
gle variable 



with 



a = sm 9, 



tan( 



(2.5) 



(2.6) 



Superexchange contributions to H can be expressed in 
the form 



Si • Sj 



{ij)\h 

+r2 
2 



A[j'+;r(n27+nj7)-l 



A 



(7) 



(2.8) 



where one recognizes a structure similar to that for su- 
perexchange in cubic vanadates f^ii^ with separation into 
a spin projection operator on the triplet state, (Si-Sj + j), 

and an operator {Si ■ Sj — j) which is finite only for 
low-spin excitations. These operators are accompanied 
by coefficients (r'i,r2,r3) which depend on the Hund ex- 
change parameter p.2p , and are given from the multiplet 
structure of d^ ions"*^ by 



1 



1 



1 



r2 = 



1 — rj 



l + 2r] 



(2.9) 



The Coulomb and Hund exchange elements deduced from 
the spectroscopic data of Zaanen and Sawatzky^ are 
U = 4.35 eV and Jh = 0.59 eV, giving a reahstic value 
of ~ 0.136 for Ti'^+ ions. For one finds^^ U = 4.98 
eV and Jh = 0.64 eV, whence rj ~ 0.13, and the values 
for ions are expected to be very similar. Finally, 
for Co3+ ions,ii U = 6.4 eV and Jh = 0.84 eV, giving 
again rj ~ 0.13. The value rj = 0.13 therefore appears to 
be quite representative for transition-metal oxides with 
partly filled t2g orbitals, whereas somewhat larger values 
have been found for systems with active Cg orbitals due 
to a stronger Hund exchange 4 

The orbital operators Aij and Bij in Eq. (|2.8p depend 
on the bond direction 7 and involve two active orbital 
colors. 



A, 
B 



(7) 

(7) 



( + T~T~ 
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( T+T^ + TrT+ 

\ ^1 01 ' *7 J7 



1 



,(7)^(7) 



(2.10) 



+ -n) "n 



(T).(7) (2.11) 
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For illustration, in the case 7 = c {{ij) || c), the orbitals 
a and b at site i are interchanged (off-diagonal hopping) 

at site j, and the electron number operator is nl'^'' = 
n-ia + nib. The quantity n^^ in Eq. (|2.8p is the number 
operator for electrons on the site in orbitals inactive for 



hopping on bond 7, rii^ — 1 ^ 



,(7) 



or riic in this example. 



For a single bond, the orbital operators in Eq. (|2.10[) 
may be written in a very suggestive form by performing 
a local transformation in which the active orbitals are 
exchanged on one bond site, specifically \a) \h) and 



\b) 



on bond 7 



40 



Then 



31 



i (7)^(7) 



(2.12) 
(2.13) 



where the scalar product in Aij is the conventional 
expression for pseudospin-1/2 variables, and the cross 
product in Bij is defined as 



7 J7" 



(2.14) 



Equations (|2.8p and (|2.12p make it clear that for a sin- 
gle superexchange bond, the minimal energy is obtained 
either by forming an orbital singlet, in which case the op- 
timal spin state is a triplet, or by forming a spin singlet, 
in which case the preferred orbital state is a triplet; we 
refer to these bond wavefunctions respectively as (os/st) 
and (ss/ot). The two states are degenerate for 77 = 0, 
while for finite Hund exchange 



(ss/ot) 



(2.15) 
(2.16) 



and the (os/st) state is favored. This propensity for sin- 
glet formation in the a = Q limit will drive much of the 
physics to be analyzed in what follows. 

Because of the off-diagonal nature of the hopping term, 
in the original electronic basis (before the local transfor- 
mation) the orbital singlet is the state 



the absence of static dimer formation. This arises due to 
frustration, and can be shown explicitly in numerical cal- 
culations, but we will not enter into this point in more 
detail here. However, we take the liberty of retaining 
the notation of the local transformation, particularly in 
Sec. IIVI when considering dimers. Because the transfor- 
mation interchanges the definitions of FO and AO con- 
figurations, we will state clearly in each section the basis 
in which the notation is chosen. 



C. Direct Exchange 

The direct exchange part is obtained by considering 
virtual excitations of active 7 orbitals on a bond {ij) \\ 7, 
which yield 

I ^ (\ ^ 3\ ^ 1- 



(y>ll7 



Si • Sj 



r2 



X nij{l ~ Tijj) -f (1 — ni^)nj 



(2r2 + r3 



5, 



}. (2.21) 



Here there are no orbital operators, but only number 
operators which select electrons of color 7 on bonds ori- 
ented along the 7-axis. When only only one active orbital 



is occupied [ni^{l 



n 



31 



)], this electron can gain energy 



— \J from virtual hopping at 77 = 0, a number which 
has only a weak dependence on the bond spin state at 
?7 > 0. When both active orbitals are occupied (rii^rij^), 
placing the two electrons in a spin singlet yields the far 
lower bond energy — J, and thus again one may expect 
much of the discussion to follow to center on dimer-based 
states of the extended system. Again the triplet spin 
excitation corresponds to the lowest energy, [U — iJii), 
and only the lower two excitations involve spin singlets 
which could minimize the bond energy. The structure 
of these terms is the same as in the ID Cg spin-orbital 
model,''* or the case of the spinel MgTi204.^*' A simpli- 
fied model for the triangular-lattice model in this limit, 
using a lowest-order expansion in 77 for the spin but not 
for orbital interactions, was introduced in Ref. i4S . 



|Vos)-^(|aa)-|6&)) 
while the orbital triplet states are 



\i^ot+) 



= l«&> 

1 

= \ba) 



{\aa) + \bb)) , 



(2.17) 

(2.18) 
(2.19) 
(2.20) 



The locally transformed basis then gives a clear analogy 
which can be used for single bonds and dimer phases in 
combination with all of the understanding gained for the 
Heisenberg model. However, we stress here that the local 
transformation fails for systems with more than 1 bond in 



D. Mixed Exchange 



Finally, the two different types of hopping channel may 



also contribute to two-step, virtual d\d}j 
tations with one off-diagonal {t) and one diagonal {t') 
process. The occupied orbitals are changed at both sites 
(Fig. [2]), and as for the superexchange term the result- 
ing effective interaction may be expressed in terms of or- 
bital fluctuation operators. To avoid a more general but 
complicated notation, we write this term only for c-axis 
bonds, 

(»i>l|c 
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where the orbital operators are 

T- T-^alc,. (2.23) 

These definitions are selected to correspond to the t~ 
pseudospin components of both operators being \bi) for 

T;^^ and \ci) for 7;|. The form of the H^^^ and terms 
is obtained from Eq. (|2.22p by a cyclic permutation of 
the orbital indices. By inspection, this type of term is 
finite only for bonds whose sites are occupied by linear 



superpositions of different orbital colors, and creates no 
strong preference for the spin configuration at small 77. 



E. Limit of vanishing Hund exchange 

In the subsequent sections we will give extensive con- 
sideration to the model of Eq. (|2.7p at 77 = 0. In this 
special case the multiplet structure collapses (spin sin- 
glet and triplet excitations are degenerate), one finds a 
single charge excitation of energy U, and the Hamiltonian 
reduces to the form 



(7) (7) 
-n} "n)" 
4 * J 



+ 2 ("«7 ~^ 



S'i'Sj ^ \Tli^Tl-j 



(2.24) 



which depends only on the ratio of superexchange to di- 
rect exchange (0 < a < 1). The first hne of Eq. (PTM]) 
makes explicit the fact that the spin and orbital sectors 
are completely equivalent and symmetrical at a = 0, at 
least at the level of a single bond. However, we will 
show that this equivalence is broken when more bonds are 
considered, and no higher symmetry emerges because of 
the color changes involved for different bond directions, 
which change the SU(2) orbital subsector. The second 
line of Eq. (|2.24p emphasizes the importance of bond oc- 
cupation and singlet formation at a = 1 (Sec. IIC). 

In the third hne of Eq. (|2.24p . the labels 7 7^ 7 re- 
fer to the two mixed orbital operators on each bond 
[Eq. (|2.23p ]. Orbital fluctuations are the only processes 
contributing to the mixed terms in this limit, where the 
spin state of the bond has no effect. We draw the at- 
tention of the reader to the fact that for the parameter 
choice a = 0.5, an electron of any color at any site has the 
same matrix element to hop in any direction. However, 
because of the different color changes involved in these 
processes, again the spin-orbital Hamiltonian does not 
exhibit a higher symmetry at this point, a result reflected 
in the different operator structures of superexchange and 
direct exchange components. 



III. LONG-RANGE-ORDERED STATES 

In this Section we study possible ordered or partially 
ordered states for the Hamiltonian of Eq. (|2.7p . As ex- 
plained in Sec. II, the parameters of the problem are 



the ratio of the direct and superexchange interactions, 
a (12. 5p . and the strength of the Hund exchange interac- 
tion, rj (|1.2p . Regarding the latter, we will discuss briefly 
the transition to ferromagnetic (FM) spin order for in- 
creasing rj in this framework. 

The first necessary step in any analysis of such an 
interacting system is to establish the energies of differ- 
ent (magnetically and orbitally) ordered states. The 
high connectivity of the triangular-lattice system sug- 
gests that ordered states will dominate, and claims of 
more exotic ground states are justifiable only when these 
are shown to be uncompetitive. The calculations in this 
Section will be performed for static orbital and spin con- 
figurations, with the virtual processes responsible for (su- 
per)exchange as the only fluctuations. In the language 
of the discussion in Sec. I, fully ordered states gain only 
potential energy at the cost of sacrificing the kinetic (res- 
onance) energy from fluctuation processes, which we will 
show in Sees. IV and V is of crucial importance here. 



A. Possible orbital configurations 

The results to follow will be obtained by first fixing 
the orbital configuration, either on every site or on par- 
ticular bonds, and then computing the spin interaction 
and optimizing the spin state accordingly. While this is 
equivalent to the converse, the procedure is more trans- 
parent and offers more insight into the candidate phases. 
We limit the number of states to ordered phases with 
small unit cells, and the orbital states to be considered 
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(e) (f) 

FIG. 3: (Color online) Schematic representation of possible 
orbital states with a single color on each site of the triangu- 
lar lattice: (a) one-color state; (b) and (c) two inequivalent 
two-color states; (d) three-sublattice three-color state; (e) 
and (f) two inequivalent three-color states. The latter two 
configurations are degenerate with similar states where the 
lines of occupied a and b orbitals repeat rather than being 
staggered along the direction perpendicular to the lines of 
occupied c orbitals. The three-sublattice state (3d) is nonde- 
generate (d = 1), states (3a), (3b), and (3e) have degeneracy 
d = 3, and states (3c) and (3f) have degeneracy d = 6. 



are enumerated in this subsection. For clarity we adopt 
the convention of Fig.[2Ic) that horizontal (c) bonds have 
diagonal (direct exchange) hopping of c orbitals, which 
are shown in blue, and off-diagonal (superexchange) hop- 
ping processes for a and h orbitals [Fig. [U^b)] , respectively 
red and green; up-slanting (a) bonds have diagonal hop- 
ping for a orbitals and off-diagonal hopping between b 
and c orbitals; down-slanting (6) bonds have diagonal 
hopping for h orbitals and off-diagonal hopping between 
a and c orbitals. All Hamiltonians and energies are func- 
tions of a and ry, as given by Eqs. (j^ . (I^?^ . 
and l|2.22p . To minimize additional notation, they will 
be quoted in this and in the next section as functions of 
the single argument a, with implicit 77-dependence con- 
tained in the parameters (7'i,r2,r3). The orbital bond 
index 7 will also be suppressed here and in Sec. IV. 

We continue to refer to the orbital type as a "color" , 
and begin by listing symmetry-inequivalent states where 




FIG. 4: (Color online) Schematic representation of possible 
orbital configurations with superpositions of (a) two orbitals 
in a two-color state, (b) three orbitals, (c) two orbitals with 
equal net weight, and (d) and (e) two orbitals with differing 
net weights of all three orbitals. State (a) has degeneracy 
d = 3, states (b) and (c) have d = 1, and the degeneracies of 
states (d) and (e) are d = Q and d = 3. 



each site has a unique color. If the same orbital is occu- 
pied at every site [Fig. Elja)], the three states with a, 6, 
or c orbitals occupied are physically equivalent (degener- 
acy is d = 3). When lines of the same occupied orbitals 
alternate along the perpendicular direction there are two 
basic possibilities, which are shown in Figs.[3I^b) and[3]Jc). 
These two-color states differ in their numbers of active 
superexchange or direct-exchange bonds, which depend 
on how the monocolored lines are oriented relative to the 
active hopping direction(s) of the orbital color. There 
is only one three-color configuration with equal occupa- 
tions, which is shown in Fig. [Sjd). 

Turning to orbital states with unequal occupations, 
motivated by the tendency of TL to favor dimer forma- 
tion in certain limits we extend our considerations to the 
possibility of a four-site unit cell [Figs. [3]Je) and[3I^f)]. 
More elaborate three-color unit cells are not considered. 
In this case the same state is obtained when the fourth 
site is occupied by electrons whose orbital color is any of 
the other three. Again this state, which breaks rotational 
symmetry, differs depending on its orientation relative to 
the active hopping axes. 
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States involving a superposition of either two or three 
orbitals at each site can be expected to allow a signif- 
icantly greater variety of hopping processes. When ei- 
ther two or three orbital states are partially occupied 
at each site (we stress that the condition of Eq. (|2.3[) 
is always obeyed rigorously), one finds the two uniform 
states represented in Figs. Hl[a) and^Kb). These denote 
the symmetric wavefunctions \ip2) — (|0a) + \4>b))/\/^ 
and IV's) = (|0a) + \<t>b) + \(j)c))/Vi at every site, where 
|(/)^) = 7^10). The remaining states shown in Fig. [5] in- 
volve only two orbitals per site, but with all three orbitals 
partly occupied in the lattice. The average electron den- 
sity per site and per orbital is 1/3 in the state of Fig.lD^c), 
while in Figs. ID^d) andSKe) it is = ^, tIq = = i. 
The latter two states are neither unique nor (for general 
interactions) equivalent to each other, and represent two 
classes of states with respective degeneracies 3 and 6. 



B. Ordered— state energies: superexchange 

Before analyzing the different possible ordered states 
for any of the model parameters, we stress that the spin 
interactions on a given bond depend strongly on the or- 
bital occupation of that bond. We begin with the pure 
superexchange model Hg (|2.8p . meaning a = 0, for which 
the question of spin and orbital singlets was addressed 
in Sec. Ill Bl Here the spin and orbital scalar products 
{Si ■ Sj) and (T; • Tj) may take only values consistent 
with long-range order throughout the system and thus 
vary between —1/4 and -1-1/4. 

For a bond on which both electrons occupy active or- 
bitals, one has the possibility of cither FO or AO states. 
For the FO state, (f, • f/) = 1/4 = (T^ x fj) and 
(Ay) — {Bij) — 1, whence the terms of Tie can be sepa- 
rated into the physically transparent form 



1 /I 



- Jri + fij^) ) ( S", • 5j + - I = 0, 



H. 



(FO) 



(0) 



Jr2 St ■ 5, - - 



(3.1) 



H. 



(FO) 



1 



J{r3 - r2) [Si- S. 



4/ ' 



specifying a net spin interaction which, because — 0, 
must be AF if any hopping processes are to occur. In 
the AO case, (f, • fj) = -1/4 = (f; x Tj) and (Ay) = 
(By ) = 0, giving 



i Jr2 (Si -Sj-^ 



(3.2) 



H. 



(AO) 



(0) = 0, 



and the spin interaction is constant at 77 = 0, with only a 
weak FM preference emerging at finite 77. We remind the 
reader here that the designations FO and AO continue 
to be based on the conventional notation^^ obtained by a 
local transformation on one bond site, and in the basis of 
the original orbitals correspond respectively to opposite 
active orbitals and to equal active orbitals. Cases where 
only one orbital is active on a bond are by definition AO, 
but do contribute a finite spin interaction 



-ijri is. -8, + -^ 



Si ' Sj 



(3.3) 



- 0, 



which again has only a weak FM tendency at ry > 0. 
Clearly, when neither electron may hop, the bond does 
not contribute a finite energy. 

We begin with the uniform, one-color orbital state of 
Fig. EJa), meaning that all bonds are AO by the defini- 
tion of the previous paragraph. In two directions both 
electrons are active, while in the third none are. The 
energy per bond is 



(3o), 



1 



FM (0) = "3"^'^1- 



(3.4) 



and the spin configuration is FM. However, an antifer- 
romagnetic (AF) spin configuration on the square lattice 
defined by the active hopping directions has energy 



4f^(0) 



- ^2) 



(3.5) 



from which one observes that all spin states are degen- 
erate at = 0. The ordered spin state spin is then FM 
for any finite 77. We note in passing that the energy per 
bond for a square lattice would have the significantly 
lower value —\J for the same Jis convention, by which 
is meant the presence of the constants -|-| and — -j in 
Eq. (|2.8p . This result is a direct refiection of the geo- 
metrical frustration of the triangular lattice, an issue to 
which we return in Sec. IVII 

The state of Fig. [2Ib) involves one set of (alternating) 
AO lines with two active orbitals and two sets of (AO) 
lines each with one active orbital. All sets of lines favor 
FM order at finite 77, with 



4M(0) = -^Jri. 



(3.6) 



Here the square-lattice state which becomes degenerate 
at 7; = 0, with 



4f^(0) 



J(ri +r2) 



(3.7) 



is more accurately described as one with two lines of AF 
spins and one of FM spins [Fig.inja)], and will be denoted 
henceforth as AFF. 
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(a) (b) 

FIG. 5: (Color online) Spin configurations minimizing the 
total energy of the superexchange Hamiltonian TLs (a = 0) 
for given fixed patterns of orbital order: (a) AFF state for 
the orbital ordering pattern of Fig. 3(c), showing how the FM 
line is selected by the direction (here b) giving zero frustration; 
(b) 60-120° ordered spin configuration minimizing the total 
energy for the orbital ordering pattern of Fig. 3(d). 

The state of Fig. [3ljc) involves one set of FO lines with 
two active orbitals, one set of lines with one active orbital, 
one half set of AO lines with two active orbitals and one 
half set of inactive lines. The two-active FO lines will 
favor AF order, while the AO and the one-active lines 
will favor FM order only at 77 > 0, giving 

4ff(0) = -^J (9ri + llr2 + Ar^) (3.8) 

from the AFF configuration, but with 2 equivalent di- 
rections for the FM line. At 77 = the energy is again 
— i J. Both _E^^p'(0) and i?^''pp(0) can be regarded as the 
energy of an unfrustrated system, in the sense that the 
spin order enforced in any one direction by the orbital 
configuration at no time denies the system the ability to 
adopt the energy-minimizing configuration in other di- 
rections. However, at finite 77 the configurations shown 
in Figs. Efb) and ^[c) will be penalized relative to the 
uniform (AO) order of Fig. [3fa) due to the presence of 
AF bonds. 

We insert here an important observation: the orbital 
state of Fig. [3]Jc) also admits the formation of ID AF 
Heisenberg spin chains on the FO (6-axis) lines. The 
energy per bond of such a state includes constant inter- 
chain contributions which are independent of the spin 
state {{Si ■ Sj) = 0) on these bonds. Of these interchain 
bonds, 1/4 are FO with two active orbitals and 1/2 have 
one active orbital. One finds 

e[o\0) = -^Jln2 (2r2 + r3) - lj(3ri +r2), (3.9) 

which gives e[^\o) ^ -0.3977 J at 77 = 0. This energy 
is significantly lower than that of an ordered magnetic 
state, a result showing that the kinetic energy gained 
from resonance processes on the chains is far more signif- 
icant than minimal potential energy gain obtainable from 
an ordering of the magnetic moments on the active inter- 
chain bonds which are active, and thus provides strong 
evidence in favor of the hypothesis that any ordered state 



will "melt" to a quantum disordered one in this system. 
We will return to this issue below. 

For the two-color superposition [Fig. Hl^a)], one set 
of bonds always has two active orbitals, but with equal 
probability of being FO or AO, while the other two sets 
of bonds have a 1/4 probability of having two active or- 
bitals, which are FO, or a 1/2 probability of having one 
active orbital (and a 1/4 probability of having none). 
Under these circumstances, the net system Hamiltonian 
can be expressed by summing over all the possible orbital 
states, although this is not necessarily a useful exercise 
when the spin state may not be isotropic. By insert- 
ing the three most obvious ordered spin states, FM, AF 
(meaning here the AF state of the triangular lattice with 
120° bond angles and {St ■ Sj) = -|) and AFF, the can- 
didate energies are 

4'm (0) = -Ijru (3.10) 

4f^(0) = 4ff(0) = -^J(5ri + 7r2 + 2r3). 

The coincidence for the results for the AF and AFF or- 
dered states in this case is an accidental degeneracy. The 
final energy E'^pjps^ = — at 77 = shows that both 
states are compromises, and it is not possible to put all 
bonds in their optimal spin state simultaneously. This 
arises because of the presence of two-active FO compo- 
nents in all three lattice directions, and will emerge as a 
quite generic feature of superposition states, albeit not 
one without exceptions. 

In general there is no compelling reason (given by Ti. 
for any value of a) to expect that two-color superposi- 
tions of this type may be favorable. While the 120° state 
of a triangular-lattice antiferromagnet is one compromise 
within a space of SU(2) operators, this type of symmetry- 
breaking is not relevant within the orbital sector, where 
there are three colors and the two-color subsector of ac- 
tive orbitals in the a — limit changes as a function of 
the bond orientation. 

In the equally weighted three-color state [Fig. El^d)], 
all bonds are FO and it is easy to show that 1/3 of them 
(arranged as isolated triangles) have two active orbitals 
while the other 2/3 have one active orbital. The two- 
active bonds favor AF order while the one-active bonds 
have only a weak preference for FM order at finite 77. 
In this case the problem becomes frustrated and is best 
resolved by a kind of AF state on the triangular lattice 
where the strong triangles have 120° angles and alternat- 
ing triangles have spins either all pointing in or all point- 
ing out [Fig. Eljb)]; then 2/3 of the intertriangle bonds 
have 60° angles while the other 1/3 have 120° angles. 
The energy of this state is 

E<-'"'\0) = J (19ri + 17r2 -f Org), (3.11) 

and ^^^^^^(O) = — ^Jat7; = 0, a value again inferior to 
the optimal energy due to the manifest spin frustration. 
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In the state of Fig. ^e), the only AO bonds (1/6 of the 
total) contain inactive orbitals. Of the remaining bonds, 
3/6 have two active FO orbitals (in all three directions) 
and 2/6 have one active orbital. Once again the system 
is composed of strongly coupled triangles, but this time 
in a square array and with strong coupling in their basal 
direction by one set of two-active FO bonds. Possible 
competitive spin-ordered states would be AF or AFF, 
with energies 



4fH0) 
^rko) 



— J(5ri + 15r2 + 6r3), (3.12) 
96 

^^J(15ri+35r2 + 16r3). 



The lowest energy is obtained for 120° AF order, with 
the frustrated value Efp{0) — for ij = 0. 

For the state in Fig. [^f) the FO bonds (1/6) and only 
1/6 of the AO bonds have two active orbitals, while the 
other 2/3 of the bonds have one active orbital. In this 
case 



E 



(3/) 
AFF 



(0) = 



-4^(15ri + 13r2 + 2r3), (3.13) 
96 

-Y^J(36ri + 17r2 + 5r3), 



leading again to an AF spin state. At 77 = one has 
Ef^\o) = —^J, i-e. relatively weaker frustration. 

Turning now to three-color superpositions, the "uni- 
form" orbital state [Fig. IHb)] is one in which on every 
bond there is a probability 2 /9 of having two active FO 
orbitals, 2/9 for two active AO orbitals, 4/9 of one active 
orbital and 1/9 of no active orbitals. The appropriately 
weighted bond interaction strengths may be summed to 
give the net interaction, which for the three spin states 
considered results in the energies 



4^(0) = 

4fH0) = 



E 



(4b) 
AFF 



(0) = 



-i:J i5n+^^2 + r3), (3.14) 
36 

-^J(12ri + 10r2 + 2r3), 



and thus the AF state is lowest, with the value 4'f'' ~ 
— II J at 77 = 0. While this orbital configuration does not 
attain the minimal energy of —-^J, it is a close competi- 
tor: although it involves every bond, the fractional prob- 
abilities of each being in a two-active state mean that 
it cannot maximize individual bond contributions. How- 
ever, we will see in Sec. HID that state (4b) lies lowest 
over much of the phase diagram (0 < a < 1) as a result 
of the contributions from mixed terms. 

For states with unequal site occupations, in Fig. [U^c) 
one has a situation where on 1/3 of the bonds (arranged 
in separate triangles) there is a 1/4 probability of two 



active FO orbitals and a 1/2 probability of one active 
orbital, while on the remaining 2/3 of the bonds there is 
a 1/4 probability of two active AO orbitals, 1/4 of two 
active FO orbitals and 1/2 of having one active orbital. 
On computing the net energies for the three standard 
spin configurations, one obtains 



4m^(0) 



4f'(0) 

4Fko) 



'24-^"^' 

^ J (25ri + 25r2 + 8r3), (3.15) 
-^J(30ri + 25r2 + 8r3), 



where the AF state with E^^p (0) = — 1| J is the lowest at 
rj = 0. However, this state is also manifestly frustrated. 

In the unequally weighted state of Fig.lH^d), the prob- 
lem is best considered once again as lines of different bond 
types. Here 1/6 of the lines have two active orbitals (1/2 
FO and 1/2 AO), 1/6 of the hues have probability 1/4 
of two active orbitals (AO) and 1/2 of one active orbital, 
1/3 of the lines have probability 1/4 of two active FO or- 
bitals, 1/4 of two active AO orbitals and 1/2 of one active 
orbital, and the remaining 1/3 of the lines have probabil- 
ity 1/4 of two active orbitals (FO) and 1/2 of one active 
orbital. The ordered spin states yield the energies 



4^(0) = 
4f (0) = 



1 

192 
1 

144 



J(25ri-H27r2 + 6r3), (3.16) 
J(21ri + 17r2+4r3), 



whence it is again the AF state, with a small degree of un- 
relieved frustration in its energy 4'f^(0) ~ which 
lies lowest at 77 = 0. 

Finally, the state of Fig. |He) has the orbital pattern of 
Fig. [3)Jd) rotated in such a way that the number of active 
orbitals in different bond directions is changed. Now 1/3 
of the bonds have probabilities 1 /4 of two active orbitals 
(AO) and 1/2 of one active orbital, while the remaining 
2/3 have probabilities 1/4 of two active orbitals (FO), 
1/4 of two active orbitals (AO) and 1/2 of one active 
orbital. The ordered-state energies are 



4m^(0) 
4fH0) 

4Fko) 



1: 



96 
1 

36 



J(15ri + 13r2 + 2r3), (3.17) 
J(6ri -I- 5r2 -f- r^), 



of which the AFF states lies lowest at = 0, achieving 

the unfrustrated value 4'ff('^) ~ That it is pos- 

sible to obtain this energy in an orbital superposition is 
because of the absence of FO bond contributions in one 
direction, which can then be chosen to be FM. 

The results of this section and the conclusions one may 
draw from them are summarized in Subsec. IlII El below. 
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C. Ordered— state energies: direct exchange 

In the limit of only direct exchange, the analysis is 
somewhat simpler. The Hamiltonian is Tid of Eq. (|2.2ip . 
and in this case a particle on any site is active in only 
one direction, which leads to the immediate observation 
that in a static orbital configuration it is never possible to 
have, on average, active exchange processes on more than 
2/3 of the bonds. For simplicity we repeat the Hamilto- 
nian for the two cases of AO order between sites, in which 
case by definition at most one of the orbitals is active, 
and FO order between sites, which is restricted to the 
case where neighboring sites have the same orbital color 
and the correct bond orientation. We stress that in this 
subsection the definitions FO and AO are entirely con- 
ventional, as the local transformation of Sec. Ill Bl is not 
relevant at a = 1, and thus the designation FO implies 
orbitals of the same color, and AO orbitals of different 
colors. One obtains the expressions 



4 



ri Si ■ Sj + - 



+r2 Si ■ S. 



\j{2r2 + n) (s^-S,-\ 



which in the 77 = limit reduce to the forms 



?^(™)(1) ^ j(s,-Sj 



(3.18) 
(3.19) 

(3.20) 
(3.21) 



It is clear (Sec. II) that for a single bond, the most fa- 
vorable state is a spin singlet, which would contribute en- 
ergy — J, but at the possible expense of placing all of the 
neighboring bonds in suboptimal states. The very strong 
preference for such singlet bonds means that any mean- 
field study of the minimal energy is incomplete without 
the consideration of dimerized (or valence-bond) states 
(Sec. ITV)) . The analysis of this section can be considered 
as elucidating the optimal energies to be gained from 
long-ranged magnetic and orbital order on these bonds, 
where the optimal energy of any one is —\J- Also as 
noted in Sec. II, any active AO bond gains an exchange 
energy {—\J) simply because it does not prevent one of 
the two particles from performing virtual hopping pro- 
cesses, and this we term "avoided blocking" . In the limit 
of zero Hund exchange, these will give a highly degener- 
ate manifold of all possible spin states, from which FM 
states are selected at finite 77. 

We begin again with one-color state of Fig.[3fa), which 
we denote henceforth as (3a). Only one set of lattice 
bonds has finite interactions, all FO, and therefore the 
system behaves as a set of AF Heisenberg spin chains 
with energy per bond 

-ijln2 (2r2-|-r3), (3.22) 



4f1d(1) 



whence -Bafid(I) = -0.2310 J at 77 = 0. 

In state (3b), the FO lines do not correspond to active 
hopping directions. The remaining two directions then 
form an AO square lattice with 



4m^(1) 



(3.23) 



This can be called a "pure avoided-blocking" energy. 
The spins are unpolarized at 77 = 0, where all bond spin 
states are equivalent, but any finite 77 will select FM or- 
der (hence the notation) . We will see in the remainder of 
this section that i? = — | J is the optimal energy obtain- 
able by a 2D ordered state in the direct-exchange limit 
(a = 1), where the net energy is generically higher than 
at a = quite simply because there are half as many hop- 
ping channels. Thus the "melting" of such ordered states 
into quasi-lD states becomes clear from the outset, and 
can be understood due to the very low connectivity of 
the active hopping network on the triangular lattice. 

In state (3c), one of the FO lines is active, and forms 
AF Heisenberg spin chains. Electrons in the other FO 
line are active only in a cross-chain direction, where their 
bonds are AO, and gain avoided-blocking energy, whence 



Et;\l) = -^J{2\n2 + l) = 



-0.1988J (3.24) 



at ry = 0. As in the preceding subsection, the coherent 
state of each Heisenberg chain is not altered by the pres- 
ence of additional electrons from other chains executing 
virtual hopping processes onto empty orbitals of individ- 
ual sites. The spin chains remain uncorrelated and only 
quasi-long-range-ordered until a finite value of 77, where 
FM spin polarization and a long-range-ordered state are 
favored. 

In the two-color superposition (4a), 1/3 of the bonds 
are inactive, while on the other 2/3 one has probabil- 
ity 1/4 of two active electrons (FO), 1/2 of one active 
(AO) and 1/4 of two inactive electrons. In this case, one 
obtains an effective square lattice on which an AF spin 
configuration is favored by the FO processes, with 



4f^(1) = 



-^J(3ri + 7r2 + 2r3), (3.25) 



so again Ep^^' {1) — at 77 = 0. 

The uniform three-color state (3d) maximizes AO 
bonds, but 1/3 of the bonds on the lattice remain in- 
active. Thus 



4m^(1) 



Jri, 



(3.26) 



and Hund exchange will select the FM spin state. 

The three-color state (3e) has FO lines oriented in 
their active direction and will, as in state (3c), form 
Heisenberg chains linked by bonds with AO order. While 
the geometry of the interchain coupling can differ de- 
pending on the orbital alignment in the inactive chains, 
it does not create a frustrated spin configuration and 
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the net energy is e'^^\i) = -B^p^l)- The state (3f) 
has only inactive FO Unes and so gains only avoided- 
blocking energy, from 2/3 of the bonds in the system, 

whence 4'4)(1) = 4=;^) (1). 

In the uniform three-color superposition (4b), every 
bond has probability 1/9 of containing two active elec- 
trons (FO), 4/9 of one active electron and 4/9 of remain- 
ing inactive. For the three different ordered spin config- 
urations considered in Subsec. HlIBI the energies are 



1 



^AF^(l) = -;^^(5ri + 5r2+r3), (3.27) 



E 



(4b) 
AFF 



72 

(1) = -^J(6ri + 5r2+r3), 



81 



and one finds the energy E'^p{l) 



-iij for the 120" 



AF state at 77 = 0. 

The three-color state (4c) is one in which 1/3 of the 
bonds (arranged on isolated triangles) have probability 
1/4 of being in a state with two active electrons and 1/2 
of containing one active electron, while on the other 2/3 
of the bonds there is simply a 1/2 probability of one 
active orbital. The respective energies are 



4m^(1) 



4f'(1) 



-E'aff(I) 



- — Jil5n + 13r2 + 2r3), (3.28) 

-^'^(18''i + 13^2 + 2r3). 
216 



At 77 



0, the energy 4m^(1) 



'■^J is minimized 



by a 120° state on the triangles, which are also isolated 
magnetically in this limit. Finite values of rj result in 
FM interactions between the triangles, and a frustrated 
problem in the spin sector which by inspection is resolved 
in favor of a net FM configuration only at large 77 (77 > 
0.23). 

Finally, the three-color states (4d) and (4e) yield two 
possibilities in the a = 1 limit, namely where one of the 
minority colors is aligned with its active direction and 
where neither is. In the former case. 



^FM 



4f^(1) 



E 



m 

AFF 



(1) = 



Jri , 

48 ' 

-^J(25ri-f 27r2-)-6r3), (3.29) 
96 



and the lowest energy £'^pp(l) = — ijat?7 = 0is given 
by the directionally anisotropic AFF spin configuration. 
This is because 1/2 of the lines, in two of the three di- 
rections, have some AF preference from their 1/4 prob- 
ability of containing two active orbitals, while the third 
direction has no preference at 77 = 0, and in any case 
favors FM spins at 77 > 0. In the latter case, the only 



AF tendencies arise along lines in a single direction, but 
avoided-blocking energy is sufficient to exclude the pos- 
sibility of a Heisenberg chain state. Here 



4m (1) = -^^^1. 



1 



E 
E 



(4e), 



(1) 

aff(1) 



AF 

(4e 



-Y^J(15ri + 13r2 + 2r3), (3.30) 
'h'^ (6?'i + 5r2-Hr3), 



whence i?^pp(l) = — ij at 77 = 0, in fact with two de- 
generate possibilities for the orientation of the FM line. 



D. 



Ordered— state energies: a = 0.5 



To illustrate the properties of the model in the presence 
of finite direct and superexchange contributions, i.e. at 
intermediate values of a, we consider the point a = 0.5. 
As shown in Sec. II, there is no special symmetry at this 
point, because the contributions from diagonal and off- 
diagonal hopping remain intrinsically different. States 
with long-ranged orbital (and spin) order at a = 0.5 
are mostly very easy to characterize, because all virtual 
processes, of both types, allowed by the given configura- 
tion are able to contribute in full to the net energy. For 
the many of the states considered in this section, the en- 
ergetic calculation for a — 0.5 is merely an exercise in 
adding the a = and a = 1 results with equal weight. 
Exceptions occur for superposition states gaining energy 
from processes contained in iJ,„ [Eq. (|2.22p ]. and are in 
fact decisive here. Because these terms involve explicitly 
a finite density of orbitals of all three colors on the bond 
in question, with the active diagonal color represented on 
both sites, only for states (4b), (4c), and (4d), but not 
(4e) [Figs. IHb-e)], will it be necessary to consider this 
contribution. 

For state (3a), in two directions both electrons are 
active by off-diagonal hopping, while in the third both 
may hop diagonally. Diagonal hopping favors an AF spin 
configuration, while the off-diagonal hopping bonds have 
only a weak preference (by Hund exchange) for FM or- 
der. The ordered-state spin solution is then a doubly 
degenerate AFF state with energy per bond 

i?(3-)(0.5) = -^J (9ri + 7r2 + 2r3), (3.31) 



0. We remind the 



giving i;(3«)(0.5) = -ij at 77 
reader that the prefactor of the superexchange and di- 
rect exchange contributions is only half as large as in 
Subsecs. Hre] and ImC] [Eq. jSj])], so the overaU effect 
of additional hopping processes in this state is in fact 
an unfrustrated energy summation. We also comment 
that, exactly at 77 = 0, there is no obvious preference 
for any magnetic order between the diagonal-hopping 
chains. Only at unrealistically large values of 77 would 
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the system sacrifice this diagonal-hopping energy to es- 
tabhsh a square-lattice FM state. At finite 77, the one- 
color orbital state represents a compromise between com- 
peting spin states preferred by the two types of hopping 
contribution. 

State (3b) has no diagonal-hopping chains, and these 
processes therefore enforce only a weak preference for a 
FM square lattice. Because the off-diagonal hopping pro- 
cesses also favor FM order at finite r) (Subsec. [ill Bp . the 
two types of contribution cooperate and one obtains 



£;(3b)(o.5) 



1 



-Jn. 



(3.32) 



State (3c) contains one half set of diagonal-hopping 
chains, which fall along one of the directions which in 
the spin state favored by the off-diagonal hopping pro- 
cesses could be FM or AF; this degeneracy will therefore 
be broken. The other half set of chains will gain only 
avoided-blocking energy from diagonal processes, which 
will take place in the FM direction and thus cause no 
frustration even at finite 77. One obtains 

i?(3-)(0.5) = -^Ji^n + 7r2 + 2r3) , (3.33) 

and thus E^^^^O.b) = -ij at 77 = from this AFF 
configuration. The additive contributions from superex- 
change and direct exchange remove the possibility that 
Heisenberg-chain states in either of the directions fa- 
vored separately by off-diagonal (Sec. IIIB) or diagonal 
(Sec. IIIC) hopping could result in an overall lowering of 
energy. 

As in Subsec. lTlI C[ in the two-color superposition (4a) 
the diagonal hopping processes are optimized by an AFF 
spin configuration. Although this is one of the degener- 
ate states minimizing the off-diagonal Hamiltonian, the 
directions of the FM lines do not match. Insertion of the 
four possible spin states yields 



E 



4m^(0.5) = 
4fH0.5) = 

AFF(0)('^-^) = 



96 



-y^J(6ri+7r2 + r3), (3.34) 



144 



J (12ri + 14r2 + Srs), 



whence the lowest final energy is E'^p\o.5) = —-^J at 
77 = 0. As noted in the previous sections for this spin 
configuration, the optimal energy for all bonds is not 
attainable within the off-diagonal hopping sector, and 
the addition of the (small) diagonal-hopping contribu- 
tion causes little overall change. 

The equally weighted three-color state (3d) has no 
lines of diagonal-hopping bonds, and in fact these con- 
tribute only avoided-blocking energy on the bonds be- 
tween the strong triangles defined by the off-diagonal 



problem, adding to the weak propensity for FM intertri- 
angle bonds arising only from the Hund exchange. The 
diagonal processes can be taken only to alter this energy, 
and not to promote any tendency towards an alteration 
of the spin state, whose energy is then 

E^^'^^O.b) = (19ri + llrs + Sr^), (3.35) 

with £;(3'')(0.5) = -iij at 77 = 0. 

State (3e) is already frustrated in the off-diagonal sec- 
tor, and diagonal-hopping processes contribute primarily 
on otherwise inactive bonds without changing the frus- 
tration conditions. For the two candidate spin configu- 
rations 



4f^(0.5) 
4Fk0.5) 



1 

96 
1 

192 



J(5ri + 127'2+4r3), 
J(llri-)-19r2-l-8r3), (3.36) 



AF-ordered state with 



a competition won by the 120° 
e'^^p\0.5) = at 77 = 0. 

State (3f) lacks active lines of diagonal-hopping pro- 
cesses, and thus the avoided-blocking energy may be 
added simply to the results for the off-diagonal sector, 
giving 



4m^(0.5) 



^FM 



i?r/^(0.5) 
4¥k0.5) 



24 
1 

192 
5 

384 



Jri, 



J (25?'i + 19r2 + 2r3), (3.37) 
J (12?'i + 5r2 -f ra), 



or a minimum of 4¥'(0-5) = ~MJ at 77 = 0. 
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^AF V^'-J; — ~96' 

In the uniform three-color superposition (4b), on ev- 
ery bond there is a probability 4/9 of having only off- 
diagonal hopping processes, 2/9 for 2 active FO orbitals 
and 2/9 for two active AO orbitals, a probability 1/9 of 
having only diagonal hopping processes, and a probabil- 
ity 4/9 of other processes. These last include the contri- 
butions from one active diagonal or off-diagonal electron, 
and mixed processes contained in the Hamiltonian 7i„i 
(I2.22P : none of these three possibilities favors any given 
bond spin configuration other than a FM orientation at 
finite 77. The net energy contributions are 



E, 



^M^(0.5) = 



E 



1^F^^(0.5) = 



4FkO-5) = 



-^J(20ri + 18r2 + 3r3), (3.38) 

--5-J (8ri +6r2 +r3), 
54 



and thus the AF state is lowest, with E^p\o.5) ~ ^j^J 
at 77 = 0. While this energy differs from that for the AFF 
spin configuration by only j^J, its crucial property is 
that it lies below the value —\J obtained by direct sum- 
mation of the superexchange and direct-exchange contri- 
butions. 
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For this orbital configuration, all three spin states gain 
a net energy of — .-^J at = from mixed processes, 
and these are sufficient, as we shall see, to reduce the 
otherwise partially frustrated ordered-state energy to the 
global minimum for this value of a. By a small extension 
of the calculation, the energy of the 120" AF spin state 
may be deduced at = for all values of a, and is given 

by 



1 



E)^p'{a) = -^J {22-na + 8y/a{l- a)). (3.39) 

Comparison with the value obtained by direct summa- 
tion, E = — i(2— Of), reveals that state (4b) is the lowest- 
lying fully spin and orbitally ordered configuration in the 
region 0.063 < a < 0.983. That this state dominates over 
the majority of the phase diagram is a direct consequence 
of its ability to gain energy from mixed processes. 

The non-uniform three-color state (4c) also presents a 
delicate competition between spin configurations of very 
similar energies. From the preceding subsections, it is 
clear that in this case diagonal and off-diagonal processes 
favor different ground states, while there will also be a 
mixed contribution from 1/3 of the bonds. The energies 
of the three standard spin configurations are 



^M^(0.5) = 



12 ^' 



4f^(0.5) 
4FkO-5) 



384 
1 

432 



J (45ri +41r2 + lOr^), (3.40) 
J (54ri +41r2 + lOrg), 



where the AF state, obtaining E^^p{0.5) = —\J is the 
lowest at 77 = 0. 

Finally, in the three-color states (4d) and (4e), which 
are composed of lines of two-color sites, this delicate bal- 
ance between different spin configurations persists. For 
configuration (4d), an AFF state with the same orien- 
tation of the FM line (along the 6-axis) is both favored 
by diagonal hopping processes and competitive for off- 
diagonal processes. With inclusion of a small contribu- 
tion due to mixed processes, the three ordered spin states 
have energies 



E~ 



(Ad), 



51 



FM (0-5) = "^'^'■l' 
r(4<i)/ 



i?X7(0.5) = - — J (85ri + 84r2 + 18r3), (3.41) 



E 



59r2 + Ura), 



from which the AFF state minimizes the energy at 77 = 
with 4f^(0.5) = ~\J. 

For state (4e), which has no mixed contribution, the 
orientations of the FM lines in the optimal AFF states do 
not match, and it is necessary, as above, to consider both 
possibilities when performing a full comparison. These 
four ordered spin states yield the energies 



4m^(0.5) 
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Jri, 



Ef;\Q.b) = -^J(15ri + 13r2 + 2r3), 



E 



(4e) 



1 



E 



AFF(o)(0-5) = ~^^J {I8r, + lir2 + 2n), 

(4e) 



AFF( 



i)(0.5) = -^J(6ri+4r2+r3), (3.42) 



among which the AF state in fact lies lowest at 77 = 0, 
achieving the weakly frustrated value -E^f^ (0-5) — ^'^J ■ 



E. Summary 

Here we summarize the results of this section in a con- 
cise form. For the superexchange model (a = 0), a con- 
siderable number of 2D ordered orbital and spin states 
exist which return the energy —^J at -q = 0. This de- 
generacy is lifted at any finite Hund exchange in favor 
of orbital states [(3a), (3b)] permitting a fully FM spin 
alignment. Most other orbital configurations introduce a 
frustration in the spin sector at small t], while some offer 
the possibility of a change of ground-state spin configu- 
ration at finite 77, where ri exceeds the r2 and contri- 
butions and begins to favor states with more FM bonds. 

However, the value E — —\J per bond remains a 
rather poor minimum for a system as highly connected 
as the triangular lattice, even if, as in the superexchange 
limit, active hopping channels exist only in two of the 
three lattice directions for each orbital color. Indeed, 
the limitations of the available ordering (potential) en- 
ergy are clearly visible from the fact that a significantly 
lower overall energy is attained in systems which aban- 
don spin order in favor of the resonance (kinetic) en- 
ergy gains available in one lattice direction. The result 
-^iD^''(0) = —0.3977 J is the single most important ob- 
tained in this section, and in a sense obviates all of the 
considerations made here for fully ordered states, man- 
dating the full consideration of 2D magnetically and or- 
bitally disordered phases. 

In the study of ordered states, it becomes clear that 
the Hund exchange acts to favor FM spin alignments at 
high 7/. Because the "low-spin" states of minimal energy 
are in fact stabilized by quantum corrections due to AF 
spin fluctuations, the lowest energies at 77 = are never 
obtained for FM states, and therefore increasing 77 drives 
a phase transition between states of differing spin and 
orbital order. We show in Fig. [6] the transitions from 
quasi-lD AF-correlated states at low 7;, for both a = {) 
and a = 1, to FM states of fixed orbital and spin order 
(3b). The transitions occur at the values 77c (0) = 0.085 
and 77c(l) = 0.097, indicating that FM ordered states 
may well compete in the physical parameter regime. We 
note again that the energies in the superexchange limit 
are lower by approximately a factor of two compared to 
the direct-exchange limit simply because of the number 
of available hopping channels. 

We note also that there is never a situation in which the 
spin Hamiltonian becomes that of a Heisenberg model on 
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FIG. 6: (Color online) Minimum energies per bond obtained 
for orbitally ordered phases, showing a transition as a func- 
tion of Hund exchange rj from quasi-lD, AF-correlated to 
FM ordered spin states. For the superexchange Hamiltonian 
"W/ of Sec. II (a = 0), the transition is from the quasi-lD 
spin state on orbital configuration (3c) [black, dashed line 
from Eq. (|3.9|l ] to the one-color orbital state (3a) [red, solid 
line from Eq. (|3.4p ]. For the direct-exchange Hamiltonian 
Tid (a = 1), the transition is from the purely Id spin state 
on the one-color orbital state (3a) [green, dot-dashed line 
from Eq. (|3.22|l ] to the two-colour, avoided-blocking state 
(3b) [blue, dotted line from Eq. (|3.23p ]. The transitions to 
FM order as obtained from the mean-field considerations of 
this section are marked by arrows. 

a triangular lattice. This demonstrates again the inher- 
ent frustration introduced by the orbital sector. However, 
the fact that the ordered-state energy can never be low- 
ered to the value i?HAF = ~^Ji which might be expected 
for a two-active FO situation on every bond, far less the 
value — ^ J which could be achieved if it were possible to 
optimize every bond in some ordered configuration, can 
be taken as a qualitative reflection of the fact that on 
the triangular lattice the orbital degeneracy "enhances" 
rather than relieves the (geometrical) frustration of su- 
perexchange interactions (Sec. IVI|) . 

The limit of direct exchange {a — 1) is found to be 
quite different: the very strong tendency to favor spin 
singlet states, and the inherent one-dimensionality of the 
model in this limit (one active hopping direction per or- 
bital color), combine to yield no competitive states with 
long-ranged magnetic order. Their optimal energy is 
very poor because of the restricted number of hopping 
channels, and coincides with the ("avoided-blocking") 
value for the model with only AO bonds, E = —^J- 
Thus these states form part of a manifold with very 
high degeneracy. However, even at this level it is clear 
that more energy, meaning kinetic (from resonance pro- 
cesses) rather than potential, may be gained by forming 
quasi-lD Heisenberg-chain states with little or no inter- 
chain coupling and only quasi-long-ranged magnetic or- 
der. Studies of orbital configurations permitting dimer- 
ized states are clearly required (Sec. IV). Finite Hund 
exchange acts to favor ordered FM configurations, which 



will take over from chain-like states at sufficiently high 
values of 77 (Fig. [5]). 

Finally, ordered states of the mixed model show a num- 
ber of compromises. At a = 0.5, where the coefficients of 
superexchange and direct -exchange are equal, some con- 
figurations are able to return the unfrustrated sum of the 
optimal states in each sector when considered separately, 
namely —\J- However, superposition states, which are 
not optimal in either limit, can redeem enough energy 
from mixed processes to surpass this value, and in fact 
the maximally superposed configuration (4b) is found 
to minimize the energy over the bulk of the phase di- 
agram. Still, the net energy of such states remains small 
compared to expectations for a highly connected state 
with three available hopping channels per orbital color. 
Because of the directional mismatch between the diago- 
nal and off-diagonal hopping sectors, no quasi-lD states 
with only chain-like correlations are able to lower the 
ordered-state energy in the intermediate regime. 

IV. DIMER STATES 

As shown in Sec. II, the spin-orbital model on a sin- 
gle bond favors spin or orbital dimer formation in the 
superexchange limit, and spin dimer formation in the 
direct-exchange limit. The physical mechanism respon- 
sible for this behavior is, as always, the fluctuation en- 
ergy gain from the highly symmetric singlet state. On 
the basis of this result, combined with our failure to find 
any stable, energetically competitive states with long- 
ranged spin and orbital order in either limit of the model 
(Sec. lIIl| ). we proceed to examine states based on dimers. 
Given the high connectivity of the triangular lattice, 
dimer-based states are not expected a priori to be capa- 
ble of attaining lower energies than ordered ones, and 
if found to be true it would be a consequence of the 
high frustration, which as noted in Sec. I has its ori- 
gin in both the interactions and the geometry. Here we 
consider static dimer coverings of the lattice, and com- 
pute the energies they gain due to inter-singlet corre- 
lations. The tendency towards the formation of singlet 
dimer states will be supported by the numerical results in 
Sec. |Vl which will also address the question of resonant 
dimer states. 



A. Superexchange model 

Motivated by the fact that the spin and orbital sectors 
in TLs (|2.8p are not symmetrical, we proceed with a sim- 
ple decoupling of spin and orbital operators. Extensive 
research on spin-orbital models has shown that this pro- 
cedure is unlikely to capture the majority of the physical 
processes contributing to the final energy, particularly 
in the vicinity of highly symmetric points of the general 
Hamiltonian. The results to follow are therefore to be 
treated as a preliminary guide, and a basis from which 



to consider a more accurate calculation of the missing 
energetic contributions. We remind the reader that the 
notation FO and AO used in this subsection is again that 
obtained by performing a local transformation on one 
site of every dimcr. As noted in Sec. Ill Bl this procedure 
is valid for the discussion of states based on individual 
dimerized bonds, where it represents merely a notational 
convenience. For FO configurations, which in the origi- 
nal basis have different orbital colors, one might in prin- 
ciple expect that, because of the color degeneracy, there 
should be more ways to realize these without frustration 
than there are to realize AF spin configurations; however, 
because of the directional dependence of the hopping, we 
will find that this is not necessarily the case (below). 

The basic premise of the spin-orbital decoupling is that 
if the spin (orbital) degrees of freedom on a dimer bond 
form a singlet state, their expectation value {Si ■ Sj) 
({Ti ■ Tj)) on the neighboring interdimer bonds will be 
precisely zero. The optimal orbital (spin) state of the 
interdimer bond may then be deduced from the effec- 
tive bond Hamiltonian obtained by decoupling. Because 
Tis depends on the number of electrons on the sites of 
a given bond which are in active orbitals, and this num- 
ber is well defined only for the dimer bonds, the effective 
Hamiltonian will be obtained by averaging over all occu- 
pation probabilities. In contrast to the pure Heisenberg 
spin Hamiltonian, here the interdimer bonds contribute 
with finite energies, and the dimer distribution must be 
optimized. A systematic optimization will not be per- 
formed in this section, where we consider only represen- 
tative dimer coverings giving the semi-quantitative level 
of insight required as a prelude to adding dimer resonance 
processes (Sec. V). 

On the triangular lattice there are three essentially 
different types of interdimer bond, which are shown in 
Fig. [7]). For a "linear" configuration [Fig.[71[a)], the num- 
ber of electrons in active orbitals on the interdimer bond 
is two; for the 8 possible configurations where one dimer 
bond is aligned with the interdimer bond under consider- 
ation [Fig.[7{b)], the number is one on the corresponding 
site and one or zero with equal probability on the other; 
for the 14 remaining configurations where neither dimer 
bond is aligned with the interdimer bond [Fig. [TI^c)], the 
number is one or zero for both sites. The number of elec- 
trons in active orbitals is then two for type ([7^), two or 
one, each with probability 1/2, for type ^h), and two, 
one or zero with probabilities 1/4, 1/2, and 1/4 for type 
©)■ 

The effective interdimer interactions for each type of 
bond can be deduced in a manner similar to the treat- 
ment of the previous section. Considering first the situa- 
tion for a bond of type ([7^) with (os/st) dimers, setting 
{Ti -Tj) ~ yields one high-spin and two low-spin terms 
which contribute 

7?(-^^)(0) = -yr,(s..S, + l), 

ij(-^^)(0) = ^^Jr,[s...S,^\), (4.1) 
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(c) 

FIG. 7; (Color online) Types of interdimer bond differing in 
effective interaction due to dimer coordination: (a) "linear", 
(b) "semi-linear", (c) "non-linear". 

Clearly h[°^''^^'' favors FM (high-spin) interdimer spin 
configurations with coefficient i, while H^^^''^^ and 
j^ipsja) ^-p (low-spin) configurations with coeffi- 

cient I (both at ?7 = 0). Because ri exceeds r2 and 
when Hund exchange is finite, one expects a critical value 
of 77 where FM configurations will be favored. Simple al- 
gebraic manipulations using all three terms suggest that 
this value, which should be relevant for a linear chain of 
(os/st) dimers, is r^c = |. In the limit 77 0, the effective 
bond Hamiltonian simplifies to 

H^°;'^\0) = lj^S,-S,-iy (4.2) 

For a bond of type ([7^) with (ss/ot) dimers, setting 
{Si ■ Sj) = on the interdimer bond yields 

Ht'^\0) = -ljr,(T,-f, + iy (4.3) 
^(ss,7a)(o^ = -ljir,-r,)(f,xf, + \y 

Here ijj^*''''^) favors AO configurations with coefficient |, 

while fjlf''''^^'' and nj^^^''^'^'' both favor FO configurations 
with coefficient j (at 77 = 0). Over the relevant range of 
Hund exchange coupling, < 77 < 1/3, there is no change 
in sign and AO configurations are always favored. The 
effective bond Hamiltonian for 77 ^ is 

H!S''^H0) = \j(f,,-f,-iy (4.4) 

For bonds of type ([TJd), when only one electron occu- 
pies an active orbital the corresponding decoupled inter- 
dimer bond Hamiltonians are, for (os/st) dimers, 

Hi°'''\0) = ^\jn[s.-S, + l), 
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(4.5) 



The final interdimer interaction is obtained by averag- 
ing over these expressions and those (|4.ip for two active 
orbitals per bond, and takes the rather cumbersome form 



H, 



(os,7b) , 



1 



,ff (0) = Y^J(r3 + 5r2-3ri)5,.5, 

-^J(9ri+5r2+r3), (4.6) 



which reduces in the hmit 77 —s- to 



(4.7) 



For (ss/ot) dimers, the situation cannot be formulated 
analogously, because if only one electron on the bond 
is active, the orbital state of the other electron has no 
influence on the hopping process, i.e. ■ Tj is not a 
meaningful quantity. The resulting expressions lead then 
to 



H. 



(ss,7b) 
off 



1 



1 



(0) - -Ji3n-r2)T,-T,~—J{r3~r2)T,xT, 



8 ' ' ' 12 

J(9ri +5r2 -fra), 



48 



which has the 77 ^ limit 



1 



(4.8) 



(4.9) 



Finally, for a bond of type ([7t), there is no contribution 
from interdimer bond states with no electrons in active 
orbitals, so the above results [(gTH |13D and (gSl \IM] 
are already sufficient to perform the necessary averaging. 
With (os/st) dimers 



1 
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J (27ri + 13r2 + 2r3), (4.10) 



which reduces in the limit 77 — s- to 



1 



(4.11) 



while for (ss/ot) dimers, 
^feff'''^(0) = ^J(3ri-r2)f, • f,-^J{rs-r2)f,xf, 



-Y^J(27ri + 13r2 + 2r3) 



which in the 77 ^ limit gives 



Hi-''^\0) = -J^T.-T,-- 



1 



(4.12) 



(4.13) 



These results have clear implications for the nearest- 
neighbor correlations in an extended system. By inspec- 
tion, systems composed of either type of dimer would 
favor AF (spin) and AO interdimer bonds, to the extent 
allowed by frustration, and "linear" [type ([7^)] bonds 
over "semi-linear" [type ^]p)] bonds over "non-linear" 
[type bond types in Fig. [71 to the extent allowed 

by geometry. Discussion of this type of state requires in 
principle the consideration of all possible dimer coverings, 
but will be restricted here to a small number of periodic 
arrays which illustrate much of the essential physics of 
extended dimer systems within this model. 

We begin by considering the periodic covering of 
Fig. Elja), a fully linear conformation (of ground-state 
degeneracy 12) whose interdimer bond types (Table I) 
maximize the possible number of bonds of type ([7]a,). 
The counterpoint shown in Fig. [5]Jb) consists of pairs 
of dimer bonds with alternating orientations in two of 
the three lattice directions, and constitutes the simplest 
configuration minimizing (to zero) the number of type- 
d?^) interdimer bonds. The coverings in Figs, ^c) and 
(d) have the same property. These configurations exem- 
plify a quite general result, that any dimer covering in 
which there are no linear configurations [type ([7^)] of 
any pair of dimers will have 1/3 type-(I7|D) bonds, and 
thus the remaining 1/2 of the bonds must be of type 
The coverings shown in Figs. [5]Ja) and (b, c, d) 
represent the limiting cases on numbers of each type of 
bond, in that any random dimer covering will have val- 
ues between these. Indeed, it is straightforward to argue 
that, in changes of position of any set of dimers within 
a covering, the creation of any two bonds of type ^jp) 
will destroy one of type ([7^) and one of type ([Th), and 
conversely. 

Having established this effective sum rule, we turn next 
to the energies of the dimer configurations. First, for 
both types of dimer [(os/st) and (ss/ot)], all states with 
equal numbers of each bond type are degenerate, sub- 
ject to equal solutions of the frustration problem. Next, 
if frustration is neglected, it is clear from Eqs. (|4.2l4.4p . 
(|4.7l4.9p . and (|4.11l4.13p . that the AF and AO energy 
values for the three bond types (obtained by substitut- 
ing — i for Si ■ Sj and Ti ■ Tj) are respectively —5 J, — f 
and —\J, which, when taken together with the sum rule, 
suggest a very large degeneracy of dimer covering ener- 
gies. 

Returning to the question of frustration, a covering of 
minimal energy is one which both minimizes the number 



TABLE I: Occurrence probabilities for bonds of each type for 
four simple periodic dimer coverings of the triangular lattice. 



configuration 


dimer 


bond (7a) 


bond (7b) 


bond (7c) 


Fig. Ha) 
Fig.iTb) 
Fig.ETc) 
Fig.Hd) 


1 

! 
! 
! 

6 


1 

6 
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(c) (d) 

FIG. 8: (Color online) Periodic dimer coverings on the trian- 
gular lattice, each representative of a class of coverings: (a) 
linear; (b) plaquette; (c) 12-site unit cell; (d) "zig-zag". 

of FM or FO bonds, and ensures that they fall on bonds 
of type ((Tj;); both criteria are equally important. For 
the dimer covering ([S^), with maximal aligned bonds, 
it is possible by using the spin (for (os/st) dimers) or 
orbital (for (ss/ot) dimers) configuration represented by 
the arrows in Fig. [Ha) to make the number of frustrated 
(FM/FO) interdimer bonds equal to 1/6 of the total. 
Bearing in mind that the 1 /6 of bonds covered by dimers 
are also FM/FO, and that at least 1/3 of bonds on the 
triangular lattice must be frustrated for coUinear spins, 
this number is an absolute minimum. [Here we do not 
consider the possibility of non-coUinear order of the non- 
singlet degree of freedom.] Further, for this configuration 
one observes that all of the FM/FO bonds already fall on 
bonds of type ([Tt), providing an optimal case with energy 

/I 1 1 1 1 1 3 \ 
= -H6 + 6-2 + 2-4 + 6-16j 
13 

= -32^ (^-l^) 

at ry = 0. This value constitutes a basic bound which 
demonstrates that a simple, static dimer covering has 
lower energy than any long-range-ordered spin or orbital 
state discussed in Sec. |IlT] in this limit {a = 0) of the 
model. 

It remains to establish the degeneracy of the ground- 
state manifold of such coverings, and we provide only 
a qualitative discussion using further examples. If al- 
ternate four-site (dimer pair) clusters in Fig. [DJa) are 
rotated to give the covering of Fig. [5Jb), the minimal 
frustration is spoiled: by analogy with Fig. [HI it is easy 
to show that, if only 1/6 of the bonds are to be frus- 




(a) (b) 

FIG. 9: (Color online) Spin or orbital configurations (black 
arrows) within (a) linear and (b) zig-zag orbital- or spin- 
singlet dimer coverings of the triangular lattice. The number 
of frustrated interdimer bonds is reduced to 1/6 of the total, 
and all are of type ((7};). This figure emphasizes that for the 
spin-orbital model, dimer singlet formation does not exhaust 
the available degrees of freedom. 

trated, then they are of type ^}p>), and otherwise 1/3 
of the bonds are frustrated if all are to be of type ([Tj:). 
On the periodic 12-site cluster [Fig. [8jc)], one may place 
three four-site clusters in each of the possible orienta- 
tions, which as above removes all bonds of type ([7^) and 
maximizes those of type ^}p). Within this cluster it is 
possible to have only four frustrated interdimer bonds 
out of 18, while between the clusters there is again an 
arrangement of the spin or orbital arrows (c/. Fig. [5]) 
with only six FM or FO bonds out of 24, for a net total 
of 1/6 frustrated interdimer bonds, of which half are of 
type (HJj). The covering of Fig.[51Jd) represents an exten- 
sion of the procedure of enlarging unit cells and removing 
four-site plaquettes, which demonstrates that it remains 
possible in the limit of no type-(I7^) bonds to reduce frus- 
tration to 1/6 of the bonds, and to bonds of type ([T}:) 
[Fig. IHb)], whence the energy of the covering is again 
^^dim(O) = -g| (14.14p . Thus it is safe to conclude that, 
for the static-dimer problem, the ground-state manifold 
for a = consists of a significant number of degener- 
ate coverings. We do not pursue these considerations 
further because of degeneracy lifting by dimer resonance 
processes, and because the energetic difii'erences between 
static dimer configurations are likely to be dwarfed by the 
contributions from dimer resonance, the topic to which 
we turn in Sec. V. 



B. Direct exchange model 

The very strong preference for bond spin singlets (the 
factor of 4 in Eq. (|2.2ip ] suggests that dimer states will 
also be competitive in this limit, even though only 1/6 of 
the bonds may redeem an energy of — J. Following the 
considerations and terminology of the previous subsec- 
tion, we note (i) that {Si ■ Sj) = on interdimer bonds 
and (ii) that in this case, interdimer bonds have energy 
— at ry = for types ([7^) and ([TJs), and for type 
(I?!;). Because any state with a maximal number (1/6) of 
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type-(17^) bonds must have only bonds of type ^) for 
the other 2/3 [states ([5^)], such a state is manifestly less 
favorable at a = 1 than those of type ([81d)"(I8J1), where 
there are no aligned pairs of dimers. In this latter case, 
the full calculation gives 



.(1) 



1 

144 



J{9n + 19r2 + Sra) 



(4.15) 



and i?dim(l) — —jJ for 77 = 0. This energy does now ex- 
ceed that available from the formation of Heisenberg spin 
chains in one of the three lattice directions (Sec. IIII C[) . 
which gave the value i?iDAF(l) — — 0.231J. 

At the level of these calculations, the manifold of de- 
generate states with this energy is very large, and its 
counting is a problem which will not be undertaken here. 
We will show in Sec. V that, precisely in this limit, no 
dimer resonance processes occur and the static dimer cov- 
erings do already constitute a basis for the description of 
the ground state. The question of fluctuations leading to 
the selection of a particular linear combination of these 
states which is of lowest energy, i.e. of a type of order- 
by-disorder mechanism, is addressed in Ref. H^. 

At finite values of the Hund exchange, this type of 
state will come into competition with the simple avoided- 
blocking states which gain, with a FM spin state, an 
energy 



-E'fm(I) 



6 



Jri, 



(4.16) 



as 2/3 of the bonds contribute with an energy of —jJri. 
The critical value of 77 required to drive the transition 
from the low-spin dimerized state to the FM state is 
found to be 



r/c = 0.1589. 



Mixed model 



(4.17) 



Because both of the endpoints, a = and a = 1, favor 
dimerized states over states of long-ranged order, it is 
natural to expect that a dimer state will provide a lower 
energy also at a = 0.5. However, we remind the reader 
that there are no intermediate dimer bases, and caution 
there is no strong reason to expect one or other of the 
limiting dimer states to be favored close to a = 0.5. By 
inspection, the energy of an a > state can be obtained 
by direct addition of the diagonal interdimer bond con- 
tributions in an (ss/ot) or (os/st) dimer state, which is 
established by pure off-diagonal hopping, because no site 
occupancies arise which allow mixed processes. For the 
same reason, no interdimer terms impede a calculation 
of the energy of an q < 1 state by summing the off- 
diagonal interdimer bond contributions in a spin-singlet 
dimer state stabilized by purely diagonal processes. We 
will not analyze the static dimer solutions for the inter- 
mediate regime in great detail, and provide only a crude 



estimate of the a = 0.5 energy by averaging over both re- 
sults at the limits of their applicability. We will make no 
attempt here to exclude other forms of disordered state 
at a = 0.5, and return to this question in Sec. V. 

For each type of bond it is straightforward to compute 
the energy gained from interdimer hopping processes of 
the type not constituting the dimer state, and the re- 
sults are shown in Table II. The first four lines give the 
energies per bond from diagonal hopping processes oc- 
curring on the bonds of the different a = dimer states, 
and conversely for the final two lines. It is clear that the 
occupations of type ([7^) bonds preclude any hopping of 
the opposite type. For a = dimer configurations, the 
interdimer diagonal hopping on ([Tja) bonds is always of 
avoided-blocking type, while on ^) bonds a blocking 
can occur, and like the other terms is evaluated using 
{Si ■ Sj). For a = 1, off-diagonal hopping on the inter- 
dimer bonds is evaluated with {Si ■ Sj) =0 between the 
spin singlets: all processes on ([TId) bonds are those for one 
active orbital; complications arise only for ([Tt) bonds, 
where an interdimer bond between parallel dimers has 
two active AO orbitals, while one between dimers which 
are not parallel has two active FO orbitals. 

At a = 0.5, the energy of an (os/st) or (ss/ot) dimer 
state augmented by diagonal hopping processes is mini- 
mized by states ([5^) and ([HJi): the interdimer bond con- 
tributions of all coverings in Fig. [8] are equal, despite the 
different type counts, so only the a = energy is decisive. 
At ?7 = 0, 
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J, (4.18) 



(4.19) 



The energy of a spin-singlet dimer state augmented by 
off-diagonal hopping is minimal in states ^p) and, cu- 
riously, (H^): although the latter has explicitly a worse 
ground-state energy than the other states shown, the 
effect of the additional hopping is strong, not least be- 
cause all interdimer type-([5}:) bonds are between parallel 



TABLE II: Additional interdimer bond energies at a — 0.5 
due respectively to (i) diagonal hopping occurring in a state 
(designated by q = 0) stabilized by off-diagonal processes 
and (ii) off-diagonal hopping in a state (a = 1) stabilized by 
diagonal processes. 



bond 



a 



a = 0, (os/st), AF 
a = 0, (os/st), FM 
a = 0, (ss/ot), AO 
a = 0, (ss/ot), FO 
a — 1, II dimers 
1, non-|| dimers 



-3^(3ri +r2) 
-3^(3ri + r-j) 
4(3ri+r-2) 



1^1 



-l^(3ri+r2) 
-i(3ri + r2) 
-i(3ri+r2) 
-^(2r2 + r3) 
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dimers. Thus at r/ = 0, 



E. 



it) 



(0.5) = - 



^^^(0.5) 



1 

2 V24 
1 /I 



1 2 
'2-3, 
11 12 1 
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' - -48^' 



1 1 1 



2 V4 3 4 



J 



2 3 2 2 3 4, 

(4.20) 



Despite the fact that these are two completely differ- 
ent expansions, it is worth noting that the two sets of 
numbers are rather similar, which occurs because the 
significantly inferior energy of the a = 1 ground state 
is compensated by the significantly greater interdimer 
bond energies available from off-diagonal hopping pro- 
cesses. However, this result also implies that no special 
combinations of diagonal and off-diagonal dimers can be 
expected to yield additional interdimer energies beyond 
this value. 

Taking the covering ([5^) as representative of the lowest 
available energy, but bearing in mind that many other 
states lie very close to this value, an average over the two 
approaches yields 



107 

384' 



-J 



(4.21) 



at 77 = 0. This number is no longer lower than the value 
obtained in Sec. HID for fully ordered states gaining en- 
ergy from mixed processes, raising the possibility that 
non-dimer-based phases may be competitive in the in- 
termediate regime, where neither of the limiting types 
of dimer state alone is expected to be particularly suit- 
able. However, we will not investigate this question more 
systematically here, and caution that the approximations 
made both in Sec. HID and here make it difficult to draw 
a definitive conclusion. 



D. Summary 

The results of this section make it clear that static 
dimer states, while showing the same energetic trend, 
are considerably more favorable than any long-range- 
ordered states (Sec. HI) over most of the phase diagram. 
As a function of a, the dimer energy increases mono- 
tonically from — H J to —jJ, and both end-point values 
also lie below the results obtained for quasi-lD spin- 
disordered states in Sec. HI. We stress that the results of 
this section are provisional in the sense that we have not 
performed a systematic exploration of all possible dimer 
coverings, but rather have focused on a small number 
of examples illustrative of the limiting cases in terms of 
interdimer bond types. More importantly, we have con- 
sidered only static dimer coverings with effective inter- 
dimer interactions: the kinetic energy contributions due 
to dimer resonance processes for all values of a < 1 are 
missing in this type of calculation. For this reason, we 



have also refrained from investigating higher-order pro- 
cesses, which may select particular dimer states from a 
manifold of static coverings degenerate at the level of the 
current considerations. Gaining some insight into the 
magnitude and effects of resonance contributions is the 
subject of the following section. 



V. EXACT DIAGONALIZATION 
A. Clusters and correlation functions 

In this Section we present results obtained for small 
systems by full exact diagonalization (ED). Because each 
site has two spin and three orbital states, the dimension 
of the Hilbert space increases with cluster size as 6-^, 
where N is the number of sites. As a consequence, we 
focus here only on systems with = 2, 3, and 4 sites: 
all three clusters can be considered as two-, three- or 
four-site segments of an extended triangular lattice, con- 
nected with periodic boundary conditions. For the single 
bond and triangle this only alters the bond energies by 
a factor of two, a rescaling not performed here, but for 
the four-site system it is easy to see that the intercluster 
bonds ensure that the system connectivity is tetrahedral. 
We will also compare some of the single-bond and tetra- 
hedron results with those for a four-site chain. Other ac- 
cessible cluster sizes (A = 5 and 6) yield awkward shapes 
which disguise the intrinsic system properties. Indeed we 
will emphasize throughout this Section those features of 
our very small clusters which can be taken to be generic, 
and those which are shape-specific. 

Given the clear tendency to dimerization illustrated in 
Sees, mil and IIVI it is to be expected that spin correla- 
tion lengths in all regimes of a are very small. To the 
extent that the behavior of the model for any parameter 
set is driven by local physics, the cluster results should be 
highly instructive for such trends as dimer formation, rel- 
ative roles of diagonal and off-diagonal hopping, dimer 
resonance processes, lifting of degeneracies both in the 
orbital sector and between states of (os/st) and (ss/ot) 
dimers, and the importance of joint spin-orbital corre- 
lations. However, generic features of extended systems 
which cannot be accessed in small clusters are those con- 
cerning questions of high system degeneracy and subtle 
selection effects favoring specific states. 

We will compute and discuss the cluster energies, de- 
generacies, site occupations, bond hopping probabili- 
ties in diagonal and off-diagonal channels (discussed in 
Sec. VC), and the spin, orbital, and spin-orbital (four- 
operator) correlation functions. All of these quantities 
will be calculated for representative values of a and 77 
covering the full phase diagram, and each contains impor- 
tant information of direct relevance to the local physics 
properties listed in the previous paragraph. Although 
the systems we study are perforce rather small, we will 
show that one may recognize in them a number of general 
trends valid also in the thermodynamic limit. 
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We introduce here the three correlation functions, 
which for a bond {ij) oriented along axis 7 are given 
respectively by 



Si] = \ ■ Sj\n), 



(5.1) 
(5.2) 



n 

\ ^{A'^Si-Sj){f,^-fj^)\n) 

n 

^ {n\Si-Sj\n) ^ {m\fij-fj^\m), (5.3) 



where d is the degeneracy of the ground state. The 
definitions of the spin {Sij) and orbital (Tij) correla- 
tion functions are standard, and we have included ex- 
plicitly all of the quantum states {|n)} which belong to 
the ground-state manifold. The correlation function Cy 
()5.3|) contains information about spin-orbital entangle- 
ment, as defined in Sec. I: it represents the difference 
between the average over the complete spin-orbital op- 
erators and the product of the averages over the spin and 
orbital parts taken separately. It is formulated in such 
a way that — means the mean-field decoupling of 
spin and orbital operators on every bonds is exact, and 
both subsystems may be treated independently from each 
other. Such exact factorizability is found^ in the high- 
spin states at large ry; its breaking, and hence the need 
to handle coupled spin and orbital correlations in a sig- 
nificantly more sophisticated manner, is what is meant 
by "entanglement" in this context. 



B. Single bond 

We consider first a single bond oriented along the c- 
axis (Fig. [TU]) . In the superexchange limit the active or- 
bitals are a and 5, while for direct exchange only the c or- 
bitals contribute in Eq. (|2.7p . As discussed in Sec. lIV Al a 
single bond gives energy — J in the superexchange model 
{a = 0) [Fig. fTUT a)]. where the ground state has degener- 
acy d = 6 at 77 = 0, from the two triply degenerate wave 
functions (ss/ot) and (os/st). At finite rj, the latter is 
favored as it permits a greater energy gain from excita- 
tions to the lowest triplet state in the (P configuration 
[Eqs. (i2T5)) and ([216]) ]. 

Although orbital fluctuations which appear in the 
mixed exchange terms in Eq. (|2.22p may in principle 
contribute at a > 0, one finds that the wave function 
remains precisely that for a = 0, i.e. (ss/ot) degenerate 
with (os/st), all the way to a = 0.5. Thus for the param- 
eter choice specified in Sec. II, the ground-state energy 
increases to a maximum of Eq = —0.5 J here [Fig. [TUT al]. 
The degeneracy d = 6 is retained throughout the regime 
a < 0.5, and only at a = 0.5 do several additional states 
join the manifold, causing the degeneracy to increase to 
d = 15. For the entire regime a £ (0.5,1], the ground 
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FIG. 10: (Color online) Evolution of the properties of a single 
bond 7 = c as a function of a at ?7 = 0: (a) energy spectrum 
(solid lines) with degeneracies as shown; (b) spin {Sij, filled 
circles), orbital {Tij, empty circles), and spin-orbital {Cij, x) 
correlations: Sij — Tij — dj — —0.25 for a < 0.5, while 
Tij — Cij — for Q > 0.5. The ground-state energy Eq is 
— J for both the superexchange {a — 0) and direct-exchange 
{a — 1) limits, and its increase between these is a result of the 
scaling convention. The transition between the two regimes 
occurs by a level crossing at a = 0.5. For a < 0.5, the 
two types of dimer wave function [(ss/ot) and (os/st)] are 
degenerate {d = 6) for resonating orbital configurations {ab}, 
while at a > 0.5, the nondegenerate spin singlet is supported 
by occupation of c orbitals at both sites [(ss/cc)]. 



state is a static orbital configuration with c orbitals occu- 
pied at both sites to support the spin singlet, and d = 1. 
The evolution of the spectrum with a demonstrates not 
only that superexchange and direct exchange are physi- 
cally distinct, unable to contribute at the same time, but 
that the two limiting wave functions are extremely ro- 
bust, their stability quenching all mixed fluctuations for 
a single bond. In this situation it is not the ground-state 
energy but the higher first excitation energy which re- 
veals the additional quantum mechanical degrees of free- 
dom active at a = compared to a = 1 [Fig. [TUTa')]. 

The spin, orbital, and composite spin-orbital corre- 
lation functions defined in Eqs. (|5.ip - (|5.3p give more 
insight into the nature of the single-bond correlations. 
The degeneracy of wavefunctions (ss/ot) and (os/st) for 
< < 0.5 leads to equal spin and orbital correlation 
functions, as shown in Fig. ilOf b). and averaging over the 



different states gives Sij = Tij 



As a singlet for 



one quantity is matched by a triplet for the other, the two 
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sectors are strongly correlated, and indeed Cij = —j, in- 
dicating an entangled ground state. However, a consider- 
ably more detailed analysis is possible. Each of the six in- 
dividual states {\n)} within the ground manifold has the 
expectation value (n\{Si-Sj){Ti^-Tjj)\ri) = — which 
we assert is the minimum possible when the spin and 
pseudospin are the quantum numbers of only two elec- 
trons. It is clear that if the operator in Cij is evaluated 
for any one of these states alone, the result is zero. En- 
tanglement arises mathematically because of the product 
of averages in the second term of Eq. (j5.3p , and physically 
because the ground state is a resonant superposition of 
a number of degenerate states. We emphasize that the 
resulting value, = —j, is the minimum obtainable 
in this type of model, reflecting the maximum possible 
entanglement. We will show in Sec. VE that this value 
is also reproduced for the Hamiltonian of Eq. (12. 7|) on 
a linear four-site cluster, whose geometry ensures that 
the system is at the SU(4) point of the ID SU(2)(g)SU(2) 
mo del 1^ 

By contrast, for a > 0.5 those states favored by su- 
perexchange become excited, and the spin-singlet ground 
state has Sij — —j. The orbital configuration is charac- 
terized by {nicHjc) = 1, a rigid order which quenches all 
orbital fluctuations (indeed, the orbital pseudospin vari- 
ables Ti^ are zero). Thus the spin and orbital parts are 
trivially decoupled, giving Cij = 0. Finally, at the tran- 
sition point a = 0.5, averaging over all 15 degenerate 
states yields Stj = -0.15, T^j = -0.10, and Cy = -0.09. 
In summary, the very strong tendency to dimer forma- 
tion in the two limits a = and a — 1 precludes any 
contribution from mixed terms on a single bond, lead- 
ing to a very simple interpretation of the ground-state 
properties for all parameters. 

C. Triangular cluster 

We turn next to the triangle, which has one bond in 
each of the lattice directions a, b, and c. Unlike the case 
of the single bond, here the spin-orbital interactions are 
strongly frustrated, in a manner deeper than and quali- 
tatively different from the Heisenberg spin Hamiltonian. 
Not only can interactions on all three bonds not be sat- 
isfied at the same time, but also the actual form of these 
interactions changes as a function of the occupied or- 
bitals. The triangle is sufficient to prove (numerically 
and analytically) the inequivalence in general of the orig- 
inal model and the model after local transformation, for 
frustration reasons discussed in Sec. Ill Bl 

We begin with the observation that the results to follow 
are interpreted most directly in terms of resonant dimer 
states on the triangle. This fact is potentially surprising, 
given that the number of sites is odd and dimer forma- 
tion must always exclude one of them, but emphasizes 
the strong tendencies to dimer formation in all param- 
eter regimes of the model. For their interpretation we 
use a VB ansatz where it is assumed that one bond is 



occupied by an optimal dimer state, minimizing its en- 
ergy, and the final state of the system is determined by 
the contributions of the other two bonds. This ansatz is 
perforce only static, and breaks the symmetry at a crude 
level, but enables one to understand clearly the effects of 
the resonance processes captured by the numerical stud- 
ies in restoring symmetries and lowering the total energy. 

Considering first the VB ansatz for the superexchange 
model, the energy — J may be gained only on a single 
bond, in one of two ways. For the bond spin state to 
be a singlet {S = 0, (ss/ot) wave function), two differ- 
ent active orbitals are occupied at both sites in one of 
the orbital triplet states. The other two bonds lower the 
total energy when the third site has an electron of the 
third orbital color, each gaining an energy of —0.25 J due 
to the orbital interactions in Eq. (|2.8p . The energy of 
the triangle is then i5vB(0) = —0.5 J per bond, and the 
cluster has a low-spin {S = ^) ground state with degen- 
eracy d — 6 from the combination of the orbital triplet 
and the spin state of the third electron. We stress that 
the location (a, b, or c bond) of the spin singlet does not 
contribute to the degeneracy because the three VB states 
are mixed within the ground state by the contributing 
off-dimer hopping processes. The same considerations 
applied to an (os/st) dimer on one of the bonds of the 
triangle shows that there is no color and spin state of 
the third electron which allows both non-dimer bonds 
to gain the energy —0.25 J simultaneously, so the cluster 
has a higher energy of — per bond. Thus the VB 
ansatz illustrates a lifting of the degeneracy between the 
two types singlet state, the physical origin of which lies 
in the permitted off-dimer fluctuation processes, and this 
will be borne out in the calculations below. However, the 
net spin state of the cluster has little effect on the esti- 
mated energy of the (os/st) case, and its high-spin ver- 
sion (5 = 1) will be a strong candidate for the ground 
state at higher values of rj. In the direct-exchange limit 
(a = 1), the VB ansatz for spin singlets again returns an 
energy £'vb(1) = ~^Jj also because only one non-dimer 
bond can contribute. Here the off-dimer processes are re- 
stricted to the third electron, which has arbitrary color 
and spin, and cannot mix the three VB states, whence 
the degeneracy is c? = 12. 

With this framework in mind, we turn to a description 
of the numerical calculations at all values of a, beginning 
with the most important results: at a = the degeneracy 
is d = 6, and hence VB resonance is confirmed, yielding 
an energy very much lower than the static estimate, at 
Eo = —0.75 J per bond [Fig. [TlTa)]. Thus strong orbital 
dynamics and positional resonance effects operate in the 
ground-state manifold. These break the (ss/ot)/(os/st) 
symmetry, but act to restore other symmetries broken in 
the VB ansatz. At a = 1, the energy and degeneracy 
from the VB ansatz are exact, showing that the orbital 
sector is classical andintroduces no resonance effects. 

Figure lllf a) shows the complete spectrum of the tri- 
angular cluster for all ratios of superexchange to di- 
rect exchange, and in the absence of Hund coupling. 
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FIG. 11: (Color online) (a) Energy spectrum per bond for 
a triangular cluster as a function of a for 77 = 0. Ground- 
state degeneracies are as indicated, with d = 6 at a = 
and d=12ata = l. The arrows mark two transitions in the 
nature of the (low-spin) ground state, which are further char- 
acterized in panels (b) and (c). (b) Spin {Sij, filled circles), 
orbital {Tij, empty circles), and spin-orbital {Cij, x) correla- 
tion functions on the c bond, (c) Average electron densities in 
the t2g orbitals at site 1 [Figs.[2jb,c)], showing nib (solid line) 
and nia = riic (dashed). The orbital labels are shown for a c 
bond. All three panels show clearly a superexchange regime 
for a < 0.32, a direct-exchange regime for a > 0.69, and an 
intermediate regime (0.32 < a < 0.69). A full description is 
presented in the text. 



Frustration of spin-orbital interactions is manifest in 
rather dense energy spectra away from the symmetric 
points, and in a ground-state energy per bond signifi- 
cantly higher than the minimal value —J. At a = 
the spectrum is rather broad, with a significant number 
of states of relatively low degeneracy due to the strong 
fluctuations and consequent mixing of VB states in this 
regime. However, even in this case the ground state is 
well separated from the first excited state. As empha- 
sized above, the ground-state energy, Eo{0) — — 0.75J, 
is quite remarkable, demonstrating a very strong energy 
gain from dimer resonance processes. By contrast, the 



value -Eo(l) = "^^^ P^r bond found at a = 1 is exactly 
equal to that deduced from the VB ansatz, demonstrat- 
ing that this wave function is exact. Here the excited 
states have high degeneracies, mostly of orbital origin, 
and thus the spectrum shows wide gaps between these 
manifolds of states; this effect is more clearly visible in 
Fig. [T2ljc). The degeneracies shown in Fig. [TlT a) are dis- 
cussed below. In the intermediate regime, many of the 
degeneracies at the end-points are lifted, leading to a 
very dense spectrum. The two transitions at a = 0.32 
and a = 0.69 appear as clear level-crossings: the inter- 
mediate ground state is a highly excited state in both 
of the limits (a = 0, 1), reinforcing the physical picture 
of a very different type of wave function dominated by 
orbital fluctuations and, as we discuss next, with little 
overt dimer character. 

The correlation functions for any one bond of the trian- 
gle are shown in Fig. [TlT a). That Sij is constant for all a 
can be understood in the dimer ansatz by averaging over 
the three configurations with one (ss/ot) or (ss) bond 
and one 'decoupled' spin on the third site, which gives 
Sij — —■J everywhere. The orbital and spin-orbital cor- 
relation functions show a continuous evolution accompa- 
nied by discontinuous changes at two transitions, where 
the nature of the ground state is altered. The orbital 
correlation function Tij = —j2 at a = may be under- 
stood as an average over the orbital triplet and the 
two non-dimer bonds (each When a increases, this 
value is weakened by orbital fluctuations, and undergoes 
a transition at a = 0.32 to a regime where orbital fluctua- 
tions dominate, and Tij is close to zero. Above a = 0.69, 
Tij becomes positive, and approaches as a — *■ 1, 
indicating that the wave function changes to the static- 
dimer limit. While Tic vanishes on the c bond here, the 
cluster average has a finite value due to the contribution 
Tij — J from the active non-singlet bond. 

The spin-orbital correlation function Cy also marks 
clearly the three different regimes of a. When a < 0.32, 
Cij has a significant negative value [Fig.fTlTb)] whose pri- 
mary contributions are given by the four-operator com- 
ponent {{Si-Sj){TiYTjj)) . By contrast, Cij is close to zero 
in the intermediate regime, increasing again to positive 
values for a > 0.69. For all a > 0.32, dj can be shown 
to be dominated by the term —SijTij in Eq. (|5.3p . while 
the four-operator contribution is small, and vanishes as 
a ^ 1. Thus entanglement, defined as the lack of fac- 
torizability of the spin and orbital sectors, can be finite 
even for vanishing joint spin-orbital dynamics. 

Further valuable information is contained in the or- 
bital occupancies at individual sites [Fig. [TlTb)]. which 
show clearly the three different regimes. Although there 
is always on average one electron of each orbital color 
on the cluster, these are not equally distributed, as each 
site participates only in two bonds and the symmetry is 
broken. A representative site, labelled 1 in Figs. [2l^b,c)] 
has only a and c bonds, and hence the electron density 
in the b orbital is expected to differ from the other two. 
The values Ub = 2/3 and ria = He = 1/6 found in the 
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regime a < 0.32 is understood readily as following from 
a 1/3 average occupation of (ab) and (be) orbital triplet 
states on the c and a bonds, respectively, and of an (ac) 
orbital triplet state on the b bond, which ensures that the 
electron at site 1 is in orbital 6 [Fig. [2l^b)]. By contrast, 
in the regime a > 0.69, only the two static orbital con- 
figurations (cc) and (bb) on the c and b bonds contribute, 
and Ua = Tic = ^, while rif, = 0; when the system is in the 
third possible spin-singlet state, with a {bb) orbital state 
on the b bond, the third electron is either a or c. Between 
these two regimes (0.32 < a < 0.69) is an extended phase 
with equal average occupancy of all three orbitals at each 
site, a potentially surprising result given the broken site 
symmetry of the cluster. While this may be interpreted 
as a restoration of the symmetry of the orbital sector by 
strong orbital fluctuations, including those due to terms 
in Hm (I2.22P , it does not imply a higher symmetry of the 
strongly frustrated interactions at a = 0.5. 

The spectra as a function of Hund coupling 77 are shown 
in Fig. [T^] for the a = and a = 1 limits, and at a = 0.5 
to represent the intermediate regime. The lifting of de- 
generacies as a function of 77 is a generic feature. States 
of higher spin are identifiable by their stronger depen- 
dence on ?7, and in all three panels a transition is vis- 
ible from a low-spin to a high-spin state. At a = 
[Fig. fT^ a)]. the large Iow-t; gap to the next excited state 
results in the transition occurring at the rather high value 
of rjc = 0.158. This can be taken as a further indication 
of the exceptional stability of the resonance-stabilized 
ground state in the low-spin sector. The degeneracy 
d = 12 of the high-spin state is discussed below. 

The transition to the high-spin state at a = 1 also 
occurs at a high critical value, rjc = 0.169 [Fig. [TWa)]. 
due in this case quite simply to the lack of competition 
for the strong singlet states on individual bonds. Only 
in the intermediate regime, 0.32 < a < 0.69, where we 
have shown already that the orbital state is quite differ- 
ent from that in either limit [Fig. [TT] , is the transition 
to the high-spin state much more sensitive to r/. The or- 
bital fluctuations in this phase occur both in the low-spin 
and the high-spin channel, making these very similar in 
energy, and the transition occurs for a = 0.5 at only 
?7c = 0.033 [Fig. m^b)]. As expected from the a = 
limit, where fluctuations are also strong, the characteris- 
tic features of this energy spectrum are low degeneracy 
and a semicontinuous nature. The location of the high- 
spin transition as a function of a may be used to draw a 
phase diagram for the triangular cluster, which has the 
rather symmetric form shown in Fig. 1131 

Yet more information complementary to that in the 
energy spectra and correlation functions can be obtained 
by considering the average "occupation correlations" for 
a bond (ij) \\ 7, 

P = {nijUj^), (5.4) 
Q = {ni-fil - rij^)) + {{1 - ni^)nj^), (5.5) 
R = ((l-n,^)(l-n,^)). (5.6) 

These probabilities {P + Q + R — 1) reflect directly the 




FIG. 12: (Color online) Energy spectra for a triangular cluster 
as a function of Hund exchange r). Energies are quoted per 
bond, and shown for: (a) a = 0, (b) a = 0.5, and (c) a = 1. 
The arrows indicate transitions at ric from the low-spin {S = 
1/2) to the high-spin [S — 3/2) ground state. The numbers 
in all panels give degeneracies for the two lowest states for 
77 < 77c and rj > rjc, respectively. 



nature of the resonance processes contributing to the en- 
ergy of the cluster states, in that they show the relative 
importance of diagonal and off-diagonal hopping in the 
ground states, and the evolution of these contributions 
with a and rj. We do not present these quantities in de- 
tail here, but only summarize the overall picture of the 
ground state whose understanding they help elucidate. 

For this summary we return to the VB framework, 
which accounts for many of the basic properties illus- 
trated in the numerical results presented above. Con- 
sidering first the low-spin states (77 = 0), at a = the 
ground state is given by one (ss/ot) dimer resonating 
around the three bonds of the cluster; the third site has 
the third color, its hopping gives a large value of Q = 1/3 
(i? = 2/3 from the pure superexchange channel) and its 
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FIG. 13: (Color online) Phase diagram of the triangular clus- 
ter in the plane {a,r]). The spin states below and above the 
transition line ric{a) are respectively spin doublet (5* — 1/2) 
and spin quartet (5* = 3/2). 



spin an addition twofold degeneracy (d = 3 x 2 = 6); 
the orbital occupation of the (ss/ot) dimer is responsible 
for the net 1/6:1/6:2/3 occupation distribution. When 
a > the state remains essentially one with a resonat- 
ing spin singlet, large Q and dominant R, but the orbital 
triplet degeneracy is lifted to 2 -f 1 and the ground-state 
degeneracy to d = 2 x 2. All quantities, including P, Q, 
and _R, undergo discontinuous changes at a ~ 0.32, and 
in this regime there is no longer strong evidence for an 
interpretation in terms of resonating spin singlets: large 
(5 — 2/3 and the equal site occupations suggest the dom- 
inance of mixed hopping processes which are not con- 
sistent with either mechanism of singlet formation. The 
retention of fourfold degeneracy across this transition is 
largely accidental, and stems from twofold spin and or- 
bital contributions. Only for a > 0.69 is a spin-singlet 
description once again valid: here P becomes significant, 
as the resonating singlet is stabilized by diagonal hopping 
where the orbital has the bond color. The third site now 
has one of two possible colors, its hopping keeps Q large, 
and its spin yields another twofold degeneracy, as do the 
orbital states, whence the net degeneracy is d = 2 x 2 x 2. 
Only at a = 1 does the spin singlet become static, while 
the third site still has cither of the other colors, yield- 
ing the symmetric result P = 1/9, Q = R = 4/9, and 
degeneracy d = 12. 

A similar description is possible in the high-spin states 
at ?7 > 7]c- At a = the (os/st) dimer is rendered static 
by the fact that hopping to the third site is now excluded 
if it has the third color, and so instead this site takes one 
of the singlet colors, a twofold degree of freedom which, 
however, does not allow singlet motion; as a consequence 
the orbital occupation is uniform (1/3:1/3:1/3), the hop- 
ping processes include contributions in the diagonal chan- 
nel {P — 1/6, Q = 1/3, R = 1/2) and the degeneracy 
is(i = 3x4 = 12. For a > the orbital singlet may 



again resonate, but the third site retains one of the sin- 
glet colors, orbital degeneracy is broken and d = 4. Once 
again strong mixed processes dominate the intermedi- 
ate regime, in which the spin state is not an important 
determining factor. Above a — 0.69 the critical value 
rjc required to overcome spin singlet formation becomes 
large again, and the high-spin state is one where avoided- 
blocking processes (large Q) dominate, while broken or- 
bital degeneracy keeps d = 4. Finally, at a = 1 one 
obtains a pure avoided-blocking state with orbital con- 
figurations acb or cba for the sites (1,2,3) of Fig. [DJc), 
and consequent degeneracy d = 4 x 2 = 8. Thus it is 
clear that the high-77 region is also one yielding interest- 
ing orbital models with nontrivial ground states, some 
including orbital singlet states. 

D. Tetrahedral cluster 

As in the case of the triangular lattice, interpretation of 
the numerical results for the tetrahedral cluster (four-site 
plaquette of the triangular lattice) is aided by considera- 
tion of the VB ansatz in the two limits of superexchange 
and direct-exchange interactions. The tetrahedral clus- 
ter can accommodate exactly two dimers, with all inter- 
dimer bonds of type ([Tt), and may thus be expected to 
favor dimer-based states by simple geometry. However, 
because the considerations and comparisons of this sub- 
section are given only for this single cluster type, any bias 
of this sort would not invalidate the results and trends 
discussed here. 

Because of the different forms and symmetries of the 
spin and orbital sectors, there is no possibility of elemen- 
tary spin-orbital operators, or of a ground-state wave 
function which is a net singlet of a higher symmetry 
group. The state with two orbital singlets on one pair 
of bonds, two spin singlets on a second pair and pure 
interdimer bonds on the third pair does exist, but is not 
competitive: the energy cost for removing the orbital sin- 
glets from the spin state maximizing their energy is by 
no means compensated by the energy gain from having 
two spin singlet bonds in an orbital state which also does 
not maximize their energy. This result may be taken as 
a further indication for the stability of dimers only in 
the forms (os/st) or (ss/ot) in this model, and states of 
shared orbital and spin singlets are not considered further 
here. We return to this point in the following subsection, 
in the context of the four-site chain. 

We discuss only the energies of the VB wave functions 
at 77 = 0. The minimal values obtainable for {Si ■ Sj) and 
{Ti~^-Tj~^) on the interdimer bonds is —1/4, corresponding 
to the AF/AO order. Thus at a = the energy per bond 
is 

Eos/sM ^ E,,f,t{0) ^ -^J, (5.7) 

with the degeneracy of the (ss/ot) and (os/st) wave func- 
tions restored as for the single bond. In the limit of direct 




exchange, the VB wave function consists of spin singlets 
with two active orbitals of the bond. The geometry of the 
cluster precludes these orbitals from being active on any 
of the interdimer bonds, as a result of which the energy 
per bond at 77 = is 
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and the ground state has degeneracy d — 3. 

The most important results for the tetrahedron, which 
we discuss in detail in the remainder of the subsection, 
are the following. At a = 0, the exact ground state 
energy is Eq — — 0.5833J: while not as large as in the 
case of the triangle (Sec. IV CI) , the resonance energy 
contribution is very significant also for an even number 
of cluster sites. The degeneracy of the numerical ground 
state, d — 6, has its origin in only one of the (ss/ot) or 
(os/st) wave functions (below), demonstrating again that 
there is no sense in which the quantum fluctuations in the 
spin and orbital sectors are symmetrical, and that the VB 
ansatz is capturing the essence of the local physics only 
at a very crude level. At a = 1, as also for the triangular 
cluster, the numerical results confirm not only the energy 
given by the VB ansatz but every detail (degeneracies, 
occupations, correlations) of this state. 

We begin the systematic presentation of results by dis- 
cussing the energy spectra at 77 = [Fig. rHT a)]. As soon 
as the degeneracies of the superexchange limit (a = 0) 
are broken, the spectrum becomes very dense, and re- 
mains so across almost the complete phase diagram until 
a level-crossing at ac = 0.92. The ground-state energy 
for all intermediate values of a interpolates smoothly to- 
wards the transition, showing an initial decrease not ob- 
served in the triangle: for the tetrahedron, mixed hop- 
ping terms make a significant contribution, leading to an 
overall energy minimum around a = 0.15. The domi- 
nance of these terms is indicated by both the extremely 
high value of ac and the steepness of the low-a curve 
where the transition to the static VB phase is finally 
reached. 

The bond correlation functions shown in Fig. [HTb) 
illustrate the effects of corrections to the VB ansatz. 
The spin correlations always have the constant value 
Sij — — J, which is the most important indication of the 
breaking of symmetry between (ss/ot) and (os/st) sectors 
at low a: this value is an average over the spin-singlet 
result —3/4 (on two bonds) and four bonds with value 0, 
and thus it is clear that (ss/ot) dimers afford more res- 
onance energy. However, the proximity of (os/st) states 
suggests that a low value of 77^, the critical Hund cou- 
pling for the transition to the high-spin state, is to be 
expected (below). 

The orbital correlations average to zero at a = 0, a 
non-trivial result whose origin lies in the breaking of 
nine-fold degeneracy within the orbital sector, and re- 
main close to this value until the transition at ac- It is 
worth noting here that = implies a higher frustra- 
tion in the orbital sector than would be obtained in the 
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FIG. 14; (Color online) (a) Energy spectrum per bond for 
a tetrahedral cluster as a function of a for — 0. Ground- 
state degeneracies are as indicated, with d = 6 at a = 
and d — 150 at a = 1. The arrow marks a transition in the 
nature of the (low-spin) ground state, (b) Spin {Sij, filled 
circles), orbital (Tij, empty circles), and spin-orbital {Cij, x) 
correlation functions on the c bond of the tetrahedral cluster 
as functions of a for rj — 0. 



spin sector for an (os/st) state {Sij = ^1^), which is due 
to the complex direction-dependence of the orbital de- 
grees of freedom. This phase is maintained across much 
of the phase diagram, with only small changes to the cor- 
relation functions, the negative value of Tij reflecting an 
easing of orbital frustration. The lack of a phase transi- 
tion throughout the region in which mixed processes are 
also important suggests that a dimer-based schematic 
picture of the ground state remains appropriate for the 
four-site system, with only quantitative evolution as a 
function of a until ac = 0.92. At a = 1, the result 
Tij = — i is the consequence of c-orbital operators on 
the interdimer a and b bonds. 

Significant spin-orbital correlations, Cij ~ —0.1 at 
a = [Fig. [Mlb)]. are found to be due exclusively to 
the four-operator term at low a. While these negative 
contributions drop steadily through most of the regime 
a < Uc, signifying a gradual decoupling of orbitals and 
spins as the static limit {a = 1) is approached, near ac 
the negative value of is again enhanced by the con- 
tribution —SijTij due to the interdimer bonds. Thus, 
as for the triangle (Sec. IV C|) . the entanglement is finite, 
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FIG. 15: (Color online) Energy spectra for a tetrahedral clus- 
ter as a function of Hund exchange rj. Energies are quoted per 
bond, and shown for: (a) a = 0, (b) a = 0.5, and (c) a = 1. 
The arrows indicate transitions from the low-spin {S — 0) to 
the high-spin (S = 2) ground state. 



complete factorization is not possible, and a finite value 
Cij = — 5^ is found even at a = 1. We note here that on 
the tetrahedron there is little information in the orbital 
occupations, which are constant (n-y — ^) over the entire 
phase diagram, demonstrating only the symmetry of this 
cluster geometry, and are therefore not shown. 

The spectra as a function of Hund coupling rj are shown 
for the three parameter choices a = 0, 0.5, and 1 in 
Fig. 1151 Once again, the spectra become very dense away 
from 77 = 0. At a = [Fig. [TST a)] high-spin states are 
found also in the low-energy sector, as a consequence of 
the near-degeneracy of (ss/ot) and (os/st) states, and 
the high-spin transition occurs at a very low value of rjc 
[Fig. [TsT a)]. The direct-exchange limit is both qualita- 
tively and quantitatively different, because the quantum 
fluctuations and the corresponding energy gains are lim- 
ited to the spin sector, making the low-spin states con- 
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FIG. 16: (Color online) Phase diagram of the tetrahedral clus- 
ter in the plane {a, rj). As for the triangular cluster, the spin 
states below and above the line rjc{a) are respectively singlet 
{S — 0) and quintet (5 = 2), with no intermediate triplet 
phase. 



siderably more stable and giving rj^ = 0.175 [Fig. [TSf c)]. 
The spin excitation gap decreases gradually with increas- 
ing 77, but until just below 77c, for all values of a, the spin 
excitation is to 5 = 1 states. However, these triplet 
states are never the ground state in the entire regime 
of T], a single transition always occurring directly into 
an S* = 2 state. In the intermediate regime represented 
by a — 0.5, the energy spectrum is so dense that indi- 
vidual states are difflcult to follow (a more systematic 
analysis of the spectra in different subspaces of is not 
presented here). The high-spin transition occurs at the 
relatively high value rjc = 0.136, due mainly to the large 
energy gains in the low-spin sector from mixed exchange. 
Further evidence for the importance of the orbital exci- 
tations in Hm (|2.22p can be found in the broadening of 
the spectrum which leads to the occurrence of quantum 
states with weakly positive energies: for both superex- 
change and direct-exchange processes, the Hamiltonians 
are constructed as products of projection operators with 
negative coefficients, so positive energies are excluded. 

The low- to high-spin transition points at all values of 
a can be collected to give the full phase diagram of the 
tetrahedron shown in Fig. [111 As shown above, in the 
superexchange limit the high-spin state lies very close 
to the low-spin ground state, and the transition to an 
S = 2 spin quintet occurs at 77c = 0.017. We comment 
here that this high-spin state is in no sense classical or 
trivial, being based on orbital singlets which are stabi- 
lized by strong orbital fluctuations, and emphasize again 
that the high-spin sector also contains a manifold of rich 
problems in orbital physics, which we will not consider 
further here. The near degeneracy of (ss/ot) and (os/st) 
states is further lifted in the presence of the mixed terms 
in Tim, raising 77c to values on the order of 0.12 across the 
bulk of the phase diagram. For no choice of parameters 
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is a spin triplet state found at intermediate values of rj. 
The reentrant behavior close to a = 0.5 is an indication 
of the importance of mixed terms in stabilizing a low-spin 
state, the tetrahedral geometry providing one of the few 
examples we have found of anything other than a direct 
competition, and hence an interpolation, between the two 
limiting cases. The rapid upturn in the limit of a — s- 1 
reflects the anomalous stability of the static VB states 
in the direct-exchange limit. The very strong asymme- 
try of the transition line in Fig. [1^] contrasts sharply with 
the near-symmetry about a = 0.5 observed for the trian- 
gle fFig. [T^ . and shows directly the differences between 
those features of the phase diagram which are universal 
and those which are effects of even or odd cluster sizes in 
a dimer-based system. 

We close our discussion of the tetrahedral cluster with 
a brief discussion of degeneracies and summary of the 
picture provided by the VB ansatz with additional res- 
onance. For the orbital occupation correlations and de- 
generacies, we begin with the low-spin sector (77 = 0). 
At a = one has two (ss/ot) VBs resonating around the 
6 bonds of the cluster, a state characterized by P = 1/6, 
Q = 1/3, and R = 1/2; however, a mixing of the orbital 
triplet states lowers the degeneracy from 9 to 0? = 6. 
For a > the state is the same, with slow evolution of 
P < 1/6, Q > 1/2, and R > 1/3, but now mixed hop- 
ping terms break all orbital degeneracies, giving d — 1. 
Only when a > 0.92 is the ground state more accurately 
characterized as one based on spin singlets of the bond 
color, with significant values of P and the restoration of 
an orbital degeneracy d = 2. As a — > 1, the diagonal hop- 
ping component is strengthened (P — > 1/3) as the pair of 
bond-colored spin singlets resonates, until at a = 1 they 
become static and the degeneracy is d = 3. 

For the high-spin states in the regime 77 > 77c, at a = 
one has two resonating (os/st) VBs, with the hopping 
channels unchanged and only the spin degeneracy d — 5. 
This state is not altered qualitatively for any a < 0.92, 
a transition value independent of rj. For 0.92 < a < 
1, orbital correlations are strongly suppressed and the 
state is characterized by hopping processes largely of the 
avoided-blocking type (one active orbital, Q dominant), 
still with d = 5. Finally, a = 1 represents the limit of a 
pure avoided-blocking state (P = 0, Q — 2/3, R — 1/3), 
where the degeneracy jumps to 150, a number which can 
be understood as 5 (spin degeneracy) x [6 (number of 
two-color states with no bonds requiring spin singlets) -I- 
24 (number of three-color states with no bonds requiring 
spin singlets)]. 



E. Four— site chain 

As a fourth and final case, we present results from 
a linear four-site cluster. While not directly relevant 
to the study of the triangular lattice, this system offers 
further valuable insight into the intrinsic physics of the 
spin-orbital model. The cluster is oriented along the c- 
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FIG. 17: (Color online) Evolution of the properties of the 
four-site chain as a function of a at = 0: (a) energy spec- 
trum and (b) spin {Sij, filled circles), orbital (Tij, empty cir- 
cles), and spin-orbital {Cij, x) correlation functions. Both 
panels show a transition occurring at a level crossing at 
a = 4/7. In panel (a), the labels show a nondegenerate 
ground state (d — 1) in both regimes, which has predimi- 
nantly spin singlet character at a > 0.571, but both spin 
and orbital singlet components at a < 0.571. In panel (b), 
Sij = Tij — Cij — — J for Q < 0.571 due to a resonating (ab) 
orbital configuration, while Tij = dj = for a > 0.571 as a 
consequence of the static c orbital configuration. 



axis with periodic boundary conditions. As for the single 
bond (Sec. VB), only the a and b orbitals contribute at 
a = 0, where indeed one finds average electron densities 
per site riia — no, — i, and ti^c = 0. Likewise, at a = 1 
only the c orbitals are occupied, with n^c = 1, a result 
dictated by the spin singlet correlations, which are fully 
developed only for complete orbital occupation. 

The energy per bond for the four-site chain in the 
superexchange limit is again — J, as for a single bond 
[Fig. [TTTa')]: somewhat surprisingly, the bonds do not 
"disturb" each other, and joint spin-orbital fluctuations 
extend over the entire chain. However, in contrast to 
a single bond, this behavior is due to only one quan- 
tum state, the SU(4) singlet. In this geometry, only 
one SU(2) orbital subsector is selected, and the result- 
ing SU(2)(X)SU(2) system is located precisely at the SU(4) 
point of the Hamiltonian.— Thus, exactly as in the SU(4) 
chain, all spin, orbital and spin-orbital correlation func- 
tions are equal, Sij = Tij = Cij = —0.25, as shown 
in Fig. I17r b) . For Sij and Py , this result may be under- 
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stood as an average over equal probabilities of singlet and 
triplet states on each bond. In more detail, the condition 
set on the correlation functions by SU(4) symmetryi^ is 
^{{Si- Sj){Tic-Tjc)) = Sij = Tij, an equality also obeyed 
by the single bond (Sec. VB). The product of Sij and 
Tij in its definition ensures the identity for Cy . The 
unique ground state is nevertheless a linear superposi- 
tion of states expressed in the spin and orbital bases, 
and has not only finite but maximal entanglement. This 
state persists, with a perfectly linear a-dependence, all 
the way to a = 1, but ceases to be the ground state 
at a = I [Fig. [TTT a)]. where there is a level-crossing 
with the a — 1 ground state (also perfectly linear). This 
latter state has a completely different, fluctuation-free 
orbital configuration, with pure c-orbital occupation at 
every site, and gains energy solely in the direct-exchange 
channel. The spins and orbitals are decoupled, Tij and 
Cij vanish, and the spin state has Sij — —0.50: this re- 
sult can be understood as an equal average over bond 
states with Si ■ Sj — ~j and —j, and matches that ob- 
tained for the four-site AF Heisenberg model with a res- 
onating VB (RVB) ground state.^ The energy at a = 1, 
Eo — —0.75 J [Fig. fTTTal]. is given directly by including 
the constant term, —jJ per bond, in the definition of the 
Hamiltonian 

The results for the linear four-site cluster demonstrate 
again the competition between superexchange and direct 
exchange. The orbital fluctuations arising due to the 
mixed exchange term, Hm (|2.22p . are responsible for re- 
moving the high degeneracies of the eigenenergies in the 
limits a = and a — 1 [Fig. fTW a)]. In fact the spectrum 
of the excited states is quasi-continuous in the regime 
around a = 0.5, but has a finite spin and orbital gap ev- 
erywhere other than the quantum critical point at a = | . 

These chain results raise a further possibility for the 
spontaneous formation at a = of a ID state not dis- 
cussed in Sec. III. A set of (for example) c-axis chains, 
with only a and b orbitals occupied in the pseudospin sec- 
tor, would create exactly the ID SU(4) model, and would 
therefore redeem an energy E = —jj per bond from the 
formation of linear, four-site spin-orbital singlets. The 
energy of the triangular lattice would receive a further, 
constant contribution from the cross-chain bonds, which 
was calculated in Eq. (j3.9|) for general 77, and hence would 
be given at ry = by 
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This energy represents a new minimum compared with 
all of the results in Sec. III. That it was obtained from 
a melting of both spin and orbital order confirms the 
conclusion that ordered phases are inherently unstable 
in this class of model, being unable to provide sufficient 
energy to compete with the kinetic energy gains avail- 
able through resonance processes. That its value is now 
lower than that obtained for a static, 2D dimer covering 
(Sec. IV) is not of any quantitative significance, given 
the results of Sec. V confirming the importance of the 
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FIG. 18; (Color online) Ground-state energy per bond as a 
function of a, obtained with r) = for a triangular cluster with 
3 bonds (blue, dashed line), and a tetrahedral cluster with 6 
bonds (red, solid line). For comparison, the energies obtained 
from the VB ansatz in the limiting cases a = and a = 
1 are shown for the triangular cluster (blue, diamonds) and 
tetrahedral cluster (red, yellow-filled, open circles); at a = 
both VB energies are the same, while at a = 1 they match 
the exact solutions. Green, upward-pointing triangles show 
the static-dimer results of Sec. IV for the extended system, 
and the black, dot-dashed line the lowest energy per bond 
obtained for fully spin and orbitally ordered phases in Sec. III. 
The violet, downward-pointing triangle shows the energy of 
the orbitally ordered but spin-disordered Heisenberg-chain 
state at a = [Eq. p.9|) ] and the open, yellow-filled square 
that of the analogous state at a = 1 [Eq. (|3.22|l ]. while the 
cross shows the energy of the spin- and orbitally disordered, 
SU(4)-chain state [Eq. (HH)]. 



positional resonance of dimers. 



F. Summary 

To summarize, we have shown in this section the re- 
sults of exact numerical diagonalization calculations per- 
formed on small clusters. Detailed analysis of ground- 
state energies, degeneracies, site occupancies and a num- 
ber of correlation functions can be used to extract valu- 
able information about the local physics of the model 
across the full regime of parameters. Essentially all of 
the quantities considered show strong local correlations 
and the dominance of quantum fluctuations of the short- 
est range, with ready explanations in terms of resonating 
dimer states. 

We draw particular attention to the extremely low 
ground-state energy of the triangular cluster, which 
shows large gains from dimer resonance. The tetrahedral 
cluster also has a very significant resonance contribution, 
although more of its ground-state energy is captured at 
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the level of a static dimer model. Such a VB ansatz 
provides the essential framework for the understanding 
of all the results obtained, even for systems with odd 
site numbers. The energies and their evolution with a 
contain some quantitative contrasts between even- and 
odd-site systems, allowing further insight concerning the 
range over which the qualitative features of the cluster 
results extend. 

Focusing in detail upon these energies. Fig. [TH] sum- 
marizes the exact diagonalization results at zero Hund 
coupling, and provides a comparison not only with the 
VB ansatz, but with all of the other results obtained 
in Sees. III-V. From bottom to top are shown: the ex- 
act cluster energies including all physical processes; the 
cluster VB ansatz, showing the importance of dimer reso- 
nance energy; the static VB ansatz for extended systems, 
suggesting by comparison with clusters the effects of reso- 
nance; the energies of "melted" states with ID spin (and 
orbital) correlations; the optimal energy of states with 
full, long-ranged spin and orbital order. 

Returning to the cluster results, their degeneracies can 
be understood precisely, and demonstrate the restora- 
tion of various symmetries due to resonance processes. 
We provide a complete explanation for all the correla- 
tion functions computed, and use these to quantify the 
entanglement as a function of a, 77 and the system size. 
There is a high-spin transition as a function of 77 for all 
values of a, which sets the basic phase diagram and es- 
tablishes a new set of disentangled orbital models at high 
V- 

The extrapolation of the cluster results to states of 
extended systems, some approximations for which are 
shown in Fig. 1181 is not straightforward, and cannot be 
expected to include any information relevant to subtle 
selection effects within highly degenerate manifolds of 
states. However, with the exception of the static-dimer 
regime around a = 1, our calculations suggest that noth- 
ing subtle is happening in this model over the bulk of 
the phase diagram, where the physics is driven by large 
energetic contributions from strong, local resonance pro- 
cesses. 



VI. RHOMBIC, HONEYCOMB, AND KAGOME 
LATTICES 

In Sec. I we alluded to the question of different sources 
of frustration in complex systems such as the spin-orbital 
model of Eq. (|2.7p . More specifically, this refers to the 
relative effects of pure geometrical frustration, as under- 
stood for AF spin interactions, and of interaction frustra- 
tion of the type which can arise in spin-orbital models 
even on bipartite lattices.'' Because the interaction frus- 
tration depends in a complex manner on system geom- 
etry, no simple separation of these contributions exists. 
In this section we alter the lattice geometry to obtain 
some qualitative results with a bearing on this separation, 
by considering the same spin-orbital model on the three 




(c) 

FIG. 19: (Color online) (a) Rhombic lattice, showing a two- 
color orbitally ordered state, (b) Honeycomb lattice, showing 
a one-color orbitally ordered state, (c) Kagome lattice, show- 
ing a three-color orbitally ordered state. 



simple lattice geometries which can be obtained from the 
triangular lattice by the removal of active bonds or sites. 

The geometries we discuss are rhombic, obtained by 
removing all bonds in one of the three triangular lat- 
tice directions [Fig. fTOT a)]. honeycomb, or hexagonal, ob- 
tained by removing every third lattice site [Fig. flgT b)]. 
and kagome, obtained by removing every fourth lattice 
site in a 2x2 pattern [Fig. [TWc)]. Simple geometrical 
frustration is removed in the rhombic and honeycomb 
cases, but for Heisenberg spin interactions the kagome 
geometry is generally recognized (from the ground-state 
degeneracy of both classical and quantum problems) to 
be even more frustrated than the triangular lattice. We 
consider only the a = and a = I limits of the model, 
and 7] = 0. We discuss the results for long-range-ordered 
states (Sec. IIIip and for static dimer states (Sec. ITVl for 
all three lattice geometries. Here we do not enter into nu- 
merical calculations on small clusters, and comment only 
on those systems for which exact diagonalization may be 
expected to yield valuable information not accessible by 
analytical considerations. 



A. Rhombic lattice 

While the connectivity of this geometry is precisely 
that of the square lattice, we refer to it here as rhombic 
to emphasize the importance of the bond angles of the 
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FIG. 20: (Color online) Rhombic lattice with (a) columnar 
and (b) plaquette dimer coverings. 



chemical structure in maintaining the degeneracy of the 
orbitals and in determining the nature of the exchange 
interactions. It is worth noting that the spin-orbital 
model (|2.7[) on this lattice may be realized in Sr2V04 
(below). In the absence of geometrical frustration, the 
spin problem created by imposing any fixed orbital con- 
figuration selected from Sec. |TTT] (Figs. [3] and [3]) is gen- 
erally rather easy to solve. Further, at jy = both FM 
and AF, and by extension AFF, spin states have equal 
energies, leading to a high spin degeneracy. 

Following Sec. IIIIl the a = energies for the majority 
of the orbitally ordered states of Fig. [3] are 



(6.1) 



per bond at = for a number of possible spin con- 
figurations, whose degeneracy is lifted (in favor of FM 
lines or planes) at finite rj. Indeed, the only exceptions 
to this rule occur for the three-color state [Fig. [31[d)] 
and for orientations of the other states which preclude 
hopping in one of the two lattice directions, whose tri- 
angular symmetry properties are broken by the missing 
bond. As noted in Sec. IIIIl for superpositions is it the 
exception rather than the rule for all hopping processes 
to be maximized, but on the rhombic lattice this is pos- 
sible for the states in Fig. ID^a) and some orientations of 
those in Figs.SJd) andSJe). 

For a = 1, the energy limit even on the triangular lat- 
tice was set rather by the number of active bonds than by 
the problem of minimizing their frustration. Similar to 
the a = case, all states where the active hopping direc- 
tion is one of the two lattice directions, plus in this case 
state (3d), can redeem the maximum energy available, 
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at 7y = 0, which is simply the avoided-blocking energy, 
for a large number of possible spin configurations. Finite 
Hund exchange favors FM spin states. 

Turning to dimerized states, the calculation of the en- 
ergy of any given dimer covering proceeds as in Sec. llVi 
namely by counting for each the respective numbers of 
bonds of types dTJi), (Hb), and dZt) [Fig. [7]. For the 
rhombic lattice, lack of geometrical frustration means 



that all interdimer bonds can be chosen to be AF/AO. 
The two most regular dimer coverings of the rhombic lat- 
tice with small unit cells may be designated as "colum- 
nar" [Fig. Unja)] and "plaquette" [Fig. dOIb)] . In both 
cases, 1/4 of the bonds are the dimers, and by inspection 
1 /4 of the interdimer bonds in the columnar state are of 
type (IZ^), while the remainder are (IZt); by contrast, the 
plaquette state has no type-(17^) bonds, 1/2 type-(I7)3) 
bonds, and the remainder are of type ^^). For a = 0, 
the energies are 
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at 77 = 0, both for (ss/ot) and for (os/st) dimers. The de- 
generacy of these two limiting cases, in the sense of maxi- 
mal and minimal numbers of types-d?^) and -(03) bonds, 
suggests a degeneracy of all dimer coverings at this level 
of analytical sophistication. Further, all of these dimer 
coverings are degenerate with all of the unfrustrated or- 
dered states at 77 = 0. The selection of a true ground 
state from this large manifold of static states (order- 
by~disorder) would hinge on higher-order processes, but 
these considerations are likely to be rendered irrelevant 
by dimer resonance (Sec. V). 

For the spin-singlet dimer states at a = 1 one finds 
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at 77 = 0, and thus that, as for the triangular lattice, the 
energy is minimized by dimer configurations excluding 
linear interdimer bonds. This remains a large manifold 
of dimer coverings, whose energy is manifestly lower than 
any of the possible orbitally ordered states in this limit 
of the model, and within which order-by-disorder is ex- 
pected to operate (Sec. V)4^ 

The considerations of this subsection, extended to fi- 
nite values of 77, may be relevant in the understanding of 
experimental results for Sr2V04. These suggest weak FM 
order, accompanied by an AO order^^ which could be 
interpreted as arising from the formation of dimer pairs. 
When the oxygen octahedra distort, the threefold degen- 
eracy of the t2g orbitals is lifted, to give a model con- 
taining only two degenerate orbitals, dyz and d^z- This 
leads to a situation with Ising-like superexchange inter- 
actions and quasi-lD hole propagation in an effective t-J 
model.— 



B. Honeycomb lattice 

The situation for the honeycomb lattice is very simi- 
lar to that for the rhombic case. Again the absence of 
geometrical frustration makes it possible to obtain the 
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FIG. 21: (Color online) Honeycomb lattice with (a) columnar 
and (b) three-way dimer coverings. 



minimal energy for a number of orbital orderings, with a 
high spin degeneracy at rj = 0. For pure superexchange 
interactions, once again 
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per bond, while in the direct-exchange limit 
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both at 77 = 0, for the same physical reasons as above. 

For dimer states, on the honeycomb lattice all inter- 
dimer bonds are by definition of type ([Zl;), and again can 
be made AF/AO because frustration is absent, so the 
energies of all dimer coverings are de facto identical. By 
way of demonstration, the two simplest regular configu- 
rations, which we label "columnar" and "three-way" , are 
shown in Fig.l^Il and, from the fact that now 1/3 of the 
bonds contain dimers, their energies are 
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per bond at a — = rj. Thus static dimer states are 
again degenerate with unfrustrated ordered states in the 



superexchange limit, and detailed consideration of kinetic 
processes would be required to deduce the lowest total 
energy. In this context, the dimer coverings shown in 
Fig. [21] exemplify two limits about which little kinetic 
energy can be gained from resonance (Fig. I2lf al. where 
large numbers of dimers must be involved in any given 
process) and in which kinetic energy gains from processes 
involving short loops [the three dimers around 2 /3 of the 
hexagons, Fig. [2lT b)] are maximized. 

At a = 1, only the dimer energy is redeemed, and this 
on 1/3 of the bonds, so 
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at 77 = for a large manifold of coverings. This energy 
is once again significantly better than any of the possi- 
ble ordered states, a result which can be ascribed to the 
low connectivity. That the ground state of the extended 
system in this limit for both the rhombic and honeycomb 
lattices involves a selection from a large number of nearly 
degenerate states suggests that numerical calculations on 
small clusters would not be helpful in resolving detailed 
questions about its nature. The same model for the hon- 
eycomb geometry in the a — I limit has been discussed 
for the S = 1 compound Li2Ru03,'^* where the authors 
invoked the lattice coupling, in the form of a structural 
dimerization driven by the formation of spin singlets, to 
select the true ground state. 



C. Kagome lattice 

The kagome lattice occupies something of a special 
place among frustrated spin systems'^ as one of the most 
highly degenerate and intractable problems in existence, 
for both classical and quantum spins, and even with 
only nearest-neighbor Heisenberg interactions. Inter- 
est in this geometry has been maintained by the dis- 
covery of a number of kagome spin systems, and has 
risen sharply with the recent synthesis of a true 5=1/2 
kagome material, ZnCu3(OII)6Cl2.^^ Preliminary local- 
probe experiments^iS show a state of no magnetic or- 
der and no apparent spin gap, whose low-energy spin 
excitations have been interpreted^ as evidence for an 
exotic spin-liquid phase. Both experimentally and the- 
oretically, kagome systems of higher spins {S — 3/2 and 
5/2) are found to have flat bands of magnetic excita- 
tions, reflecting the very high degeneracy of the spin 
sector While no kagome materials are yet known with 
both spin and orbital degrees of freedom, Maekawa and 
coworkers'**^' have considered the itinerant electron sys- 
tem on the triangular lattice for a = (actually for the 
motion of holes in Naa;Co02), demonstrating that the 
combination of orbital, hopping selection, and geometry 
leads to any one hole being excluded from every fourth 
site, and thus moving on a system of four interpenetrat- 
ing kagome lattices. 
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(a) (b) 

FIG. 22: (Color online) Kagome lattice with unequally 
weighted two-color states oriented (a) with and (b) against 
the lattice direction corresponding to the majority orbital 
color. 

Considering first the energies per bond for states of 
long-ranged spin and orbital order, in a number of cases 
the values for the kagome lattice are identical to those of 
the triangular lattice. This is easy to show by inspection 
for the one-color state (3a), and for the superposition 
states (4a), (4b), and (4c), where bonds of all types are 
removed in equal number. However, for the less symmet- 
rical orbital color configurations a more detailed analysis 
of the type performed in Sec. IIIII is required, and yields 
provocative results. The two simple possibilities for or- 
dered two-color states with a single color per site are 
shown in Fig. [521 and differ only in the orientation of 
the continuous lines (the majority color) relative to the 
active orbitals. These can be considered as the kagome- 
lattice analogs of states (3b) and (3c), as well as of (3e) 
and (3f). 

When the lines of c-orbitals are aligned with the c-axis 
[Fig. dHa)], this direction is inactive at a = 0, and only 
the other two directions contribute, one with two active 
FO orbitals, mandating an AF spin state to give energy 
— \j per bond, and the other with energy —\J and no 
strong spin preference, whence 

i?(k3b)(0) = -^J (6.9) 

at 77 = for sets of unfrustrated AF chains. By contrast, 
when the lines of c-orbitals fall along the 6-direction 
[Fig. dUb)], the a = problem contains one FO and 
one AO line each with two active orbitals, and one line 
with one active orbital. Only the first requires AF spin 
alignment, while the other two lines are not frustrating, 
with the result that an energy 

^^(k3c)(0) = -^J (6.10) 

can be obtained. This value is lower than that on the 
triangular lattice, showing that for the class of models 
under consideration, where not all hopping channels are 
active in all directions, a system of lower connectivity can 
lead to frustration relief even when its geometry remains 
purely that of connected triangles. 




(a) (b) 

FIG. 23: (Color online) Kagome lattice with two different, 
equally weighted three-color states: (a) two-color lines ori- 
ented such that only one superexchange channel, plus the di- 
rect exchange channel, is active on every bond, (b) two-color 
lines oriented such that all superexchange channels are active, 
but no direct exchange channels. 



With this result in mind, we consider again the possi- 
bilities offered by different three-color states, specifically 
those shown in Fig. [23l With reference to the superex- 
change problem, the state in Fig. [TW c). which by anal- 
ogy with (3d) we denote as (k3d), contains only a small 
number of remnant triangles and isolated bonds still with 
two active orbitals. However, the state (k3dl), shown in 
Fig. I^HTa) is that which ensures that no such bonds re- 
main, and every single bond of the lattice has one active 
superexchange channel. The state (k3d2) in Fig. I^^T b) 
is that in which every single bond of the lattice has two 
active (FO) superexchange channels: this possibility can 
be realized for the kagome geometry, at the cost of creat- 
ing a frustrated magnetic problem requiring a 120° spin 
state to minimize the energy, 

£^(k3d)(0) = -^J, (6.11) 
%3di)(0) = -^J, (6.12) 
i?(k3d2)(0) = -^J. (6.13) 

Thus one finds that lower energies than the value — 
per bond, which was the lower bound for fully (or- 
bitally and spin-)ordered states on the triangular lattice, 
are again possible for three-color ordered states. How- 
ever, the residual spin frustration means that the lowest 
ordered-state energy on the kagome lattice is given by the 
unfrustrated, two-color AFF state, i?(k3c)(0) — —^J- 

We present briefly the energies of the same states at 
a = 1 , where only a maximum of one hopping channel per 
bond can be active, and as noted above this is generally a 
stricter energetic limit than any frustration constraints. 
The results at r] = are 

%3b)(l) = -Jj (6.14) 
for an AFF state gaining most of its energy from the 
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c-axis chains, and 
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due to the dearth of active orbitals in this orientation. 
Similarly, by counting active orbitals in the three-color 
states, 
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and it is the state of Fig. I^Sf a) which achieves the un- 
frustrated value — ^ J by permitting one active hopping 
channel on every bond of the kagome lattice. 

We will not discuss the orbital superposition states 
which are the analogs of (4d) and (4e), noting only 
that these present again two different possibilities on the 
kagome lattice, depending on the orientation of the ma- 
jority lines. Even with the frustration relief offered by 
this geometry for the type of model under consideration, 
superposition states contain too many hopping channels 
for all to be satisfied simultaneously, and it is not possi- 
ble to equal the energy values found respectively for the 
configurations in Figs. ISST a) and (b) at a = 1 and a = 0. 

It remains to consider dimer states on the kagome lat- 
tice, as these have been of equal or lower energy for ev- 
ery case analyzed so far. The set of nearest-neighbor 
dimer coverings of the kagome lattice is large, and for 
the S = 1/2 Heisenberg model in this geometry the spin 
singlet manifold has been proposed as the basis for an 
RVB description.^ Two dimer coverings degenerate at 
the level of the current treatment are shown in Fig. [Ml 

Dimer coverings of the kagome lattice have the prop- 
erty that 3/4 of the triangles contain one dimer. In 
this case, the other bonds of the triangle are interdimer 
bonds, one of which is of type (7b) while the other is of 
type (7c). The other 1/4 of the triangles, known^° as 
"defect triangles" , have no dimers, and their three bonds 
are either all of type (7b), with probability 1/4, or one 
each of types (7a), (7b), and (7c), with probability 3/4. 
The frustration of the system is contained in the problem 
of minimizing the number of FM/FO interdimer bonds; 
this exercise is complex and no solution is known, so only 
an upper bound will be estimated here. 

The bonds of a defect triangle connect three differ- 
ent dimers, and so one (or all three) must be FM/FO. A 
hexagon of the kagome lattice with no dimers on its bonds 
is surrounded by six non-defective triangles, one with one 
dimer by one defective neighbor, with two dimers two, 
and a hexagon with three dimers shares its non-dimer 
bonds with three defect triangles. Hexagons with odd 
dimer numbers must create a FM/FO bond between at 
least one pair of dimers, and it is reasonable to place this 
bond on the defect triangle (s) where an energy cost is al- 
ready incurred. We note immediately that the cost of re- 
versing the type-(7a) bond, \j (Sec. IVA), exceeds that 




FIG. 24: (Color online) Kagome lattice with two different 
dimer coverings, (a) and (b). In both examples, only two of 
the twelve triangles shown explicitly on the cluster are "defec- 
tive" (contain no dimer), but the reader may notice that many 
of the next twelve triangles adjoining the boundary must also 
be so. 



of reversing both interdimer bonds of a non-defective tri- 
angle, which is -\- j^J . As a consequence, we take this 
cost, which is equal to that of reversing both a non- 
defective triangle and the weakest bond of the defect tri- 
angle, to be an upper bound on the effect of frustration. 
The net energy of a dimer state for a = = 77 is then 
estimated to be 
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This is a very large number for the kagome lattice, ex- 
ceeding even the value — \J per bond (which, however, is 
of no special significance here). Thus we find that dimer 
states in this type of model are strongly favored, gaining 
a very much higher energy than even the best ordered 
states. Qualitatively, the dimer energy shares with the 
ordered-state energy the feature that it is considerably 



36 



better than anything obtainable for the triangular lattice. 
This implies that the reduced connectivity of the lattice 
geometry for a model where the orbital degeneracy pro- 
vides a number of mutually exclusive hopping channels 
makes it easier to find states where every remaining bond 
can support a favorable hopping process without strong 
frustration. 

Applying all of the above geometrical considerations 
to the direct-exchange model {a = 1), where there is no 
frustration problem between the spin singlets, one finds 
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at ry = 0. Once again this energy is significantly lower 
than the value i?dim(l) = obtained for the triangu- 
lar lattice in Eq. (|4.15p . demonstrating that the multi- 
channel spin-orbital model of the type considered here is 
less frustrated in the kagome geometry. 

We comment in closing that the dimer energies we have 
estimated are only those of static VB configurations, and, 
away from a ~ I, the possibility remains of a signifi- 
cant resonance energy gain from quantum fluctuations 
between these states (c/. Sec. V). Numerical calculations 
on small clusters of sufficient size (here at least 6 sites for 
a unit cell) would be helpful in this frustrated case. 

To summarize this section, the spin-orbital model on 
bipartite lattices appears to present competing ordered 
and dimerized states with the prospect of high degenera- 
cies. Among "frustrated" systems (in the sense of being 
non-bipartite), the kagome lattice provides an example 
where geometrical and orbital frustration effects cancel 
partially, affording favorable dimerized solutions. Thus, 
while it is possible to ascribe some of the frustration ef- 
fects we have studied in the triangular lattice to a purely 
geometrical origin, for more complex models it is in gen- 
eral necessary to extend the concept of "geometrical frus- 
tration" beyond that applicable to pure spin systems. 



VII. DISCUSSION AND SUMMARY 

We have considered a spin-orbital model representa- 
tive of a strongly interacting electron system with 
the cubic structural symmetry of edge-sharing metal- 
oxygen octahedra, conditions which lead to a triangular 
lattice of magnetic interactions between sites with un- 
broken, threefold orbital degeneracy. We have elucidated 
the qualitative phase diagram, which turns out to be very 
rich, in the physical parameter space presented by the ra- 
tio (a) of superexchange to direct-exchange interactions 
and the Hund exchange (ry). 

Despite the strong changes in the fundamental nature 
of the model Hamiltonian as a function of a and i], a 



number of generic features persist throughout the phase 
diagram. With the exception of the ferromagnetic phases 
at high 77, which effectively suppresses quantum spin fluc- 
tuations (below) , there is no long-ranged magnetic or or- 
bital order anywhere within the entire parameter regime. 
This shows a profound degree of frustration whose origin 
lies both in the geometry and in the properties of the 
spin-orbital coupling; a qualitative evaluation of these 
respective contributions is discussed below. 

All of the phases of the model show a strong preference 
for the formation of dimers. This can be demonstrated 
in a simple, static valence-bond (VB) ansatz, and is re- 
inforced by the results of numerical calculations. The 
static ansatz is already an exact description of the direct- 
exchange limit, a = 1, and gives the best analytic frame- 
work for understanding the properties of much of the re- 
mainder of the phase diagram. The most striking single 
numerical result is the prevalence of VB states even on a 
triangular cluster, and the underlying feature reinforced 
by all of the calculations is the very large additional "ki- 
netic" contribution to the ground-state energy arising 
from the resonance of VBs due to quantum fluctuations. 
It is this resonance which drives symmetry restoration in 
some or all of the spin, orbital, and translational sectors 
over large regions of the parameter space. The sole ex- 
ception to dimerization is found at high 77 and around 
a = 1, where the only mechanism for virtual hopping is 
the adoption of orbital configurations which permit one 
orbital to be active ("avoided blocking"). 

The "most exotic" region of the phase diagram is that 
at small a and rj, and this we have assigned tentatively 
as an orbital liquid. In this regime, quantum fluctua- 
tions are at their strongest and most symmetrical, and 
every indication obtained from energetic considerations 
of extended systems, and from microscopic calculations 
of a range of local quantities on small clusters, suggests 
a highly resonant, symmetry-restored phase. While this 
orbital liquid is in all probability (again from the same 
indicators) based on resonating dimers, an issue we dis- 
cuss in full below, we cannot exclude fully the possibility 
of a type of one-dimensional physics: short, fluctuating 
segments of frustration-decoupled spin or orbital chains, 
whose character persists despite the high site coordina- 
tion. It should be stressed here that the point (a, ry) = 
(0,0) is not in any sense a parent phase for exotic states in 
the rest of the phase diagram: mixed and direct exchange 
processes are qualitatively different elements, which in- 
troduce different classes of frustrated model at flnite a. 
While the matter is somewhat semantic, we comment 
only that one cannot argue for the point a = 0.5 being 
"more exotic" than a = despite having the maximal 
number of equally weighted hopping channels, because it 
does not possess any additional symmetries which man- 
date qualitative changes to the general picture. In this 
sense, the limit a — 1 serves as a valuable fixed point 
which is understood completely, and yet is still domi- 
nated by the purely quantum mechanical concept of sin- 
glet formation. 
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One indicator which can be employed to quantify "how 
exotic" a phase may be is the entanglement of spin and 
orbital degrees of freedom. We define entanglement as 
the deviation of the spin and orbital sectors from the 
factorized limit in which their fluctuations can be treated 
separately. We compute a spin-orbital correlation func- 
tion and use it to measure entanglement, finding that 
this is significant over the whole phase diagram. Qual- 
itatively, entanglement is maximal around the superex- 
change limit, which is dominated by dimers where sin- 
glet formation forces the other sector to adopt a local 
triplet state. However, for particular clusters and dimer 
configurations, the high symmetry may allow less entan- 
gled possibilities to intervene exactly at a = 0. The 
direct-exchange limit, a = 1, provides additional in- 
sight into the entanglement definition: the four-operator 
spin-orbital correlation function vanishes, reflecting the 
clear decoupling of the two sets of degrees of freedom at 
this point, but the finite product of separate spin and 
orbital correlation functions violates the factorizability 
condition. 

This preponderance of evidence for quantum states 
based on robust, strongly resonating dimers implies fur- 
ther that the (spin and orbital) liquid phase is gapped. 
Such a state would have only short-ranged correlation 
functions. However, these gapped states are part of a 
low-energy manifold, and for the extended system we 
have shown that this consists quite generally of large 
numbers of (nearly) degenerate states. The availability 
of arbitrary dimer rearrangements at no energy cost has 
been suggested to be sufficient for the deconfinement of 
elementary S = 1/2 (and by analogy T — 1/2) excita- 
tions with fractional statistics. »1 However, the spinons 
(orbitons) are massive in such a model, in contrast to 
the properties of algebraic liquid phases j^'- 

A low-spin to high-spin transition, occurring as a func- 
tion of 7], is present for all values of a. The quantitative 
estimation of rjc in the extended system remains a prob- 
lem for a more sophisticated analysis. At the qualitative 
level, large 77 can be considered to suppress quantum spin 
fluctuations by promoting parallel-spin (ferromagnetic) 
intermediate states on the magnetic ions. However, even 
when this sector is quenched, the orbital degrees of free- 
dom remain frustrated, and contain non-trivial problems 
in orbital dynamics. In the superexchange (low-a, high- 
77) region, frustration is resolved by the formation of or- 
bital singlet (spin triplet) dimers, whose resonance min- 
imizes the ground-state energy. The frustration in the 
direct-exchange (high-a, high-77) region is resolved by 
avoided-blocking orbital configurations, and order-by- 
disorder effects are responsible for the selection of the 
true ground state from a degenerate manifold of possi- 
bilities; this is the only part of the phase diagram not 
displaying dimer physics. Thus the ferromagnetic orbital 
models in both limits exhibit a behavior quite different 
from that of systems with only S = 1/2 spin degrees of 
freedom on the triangular lattice. 

We have commented on both geometry and spin- 



orbital interactions as the origin of frustration in the 
models under consideration. However, a statement such 
as "on the triangular lattice, geometrical frustration en- 
hances interaction frustration for spin-orbital models" 
must be qualified carefully. We have obtained anecdotal 
evidence concerning such an assertion in Sec. VI by con- 
sidering other lattice geometries, and find that indeed the 
same model on an unfrustrated geometry appears capa- 
ble of supporting ordered states; however, the interplay 
of the two effects is far from direct, as the kagome lattice 
presents a case where dimer formation acts to reduce the 
net frustation. Quite generally, spin-orbital models con- 
tain in principle more channels which can be used for re- 
lieving frustration, but the exact nature of the coupling of 
spin and orbital sectors may result in the opposite effect. 
Specific data characterizing mutual frustration can be ob- 
tained from the spin and orbital correlations computed 
on small clusters: as shown in Sec. V, for the triangular 
lattice there are indeed regimes where, for example, the 
effective orbital interactions enforced by the spin sector 
make the orbital sector more frustrated (higher Ty ) than 
would be the analogous pure spin problem (measured by 
S'ij), and conversely. 

We comment briefly on other approaches which might 
be employed to obtain more insight into the states of the 
extended system, with a view to establishing more deflni- 
tively the nature and properties of the candidate orbital 
liquid phase. More advanced numerical techniques could 
be used to analyze larger unit cells, but while Lanczos 
diagonalization, contractor renormalization^ or other 
truncation schemes might afford access to systems two, 
or even four, times larger, it seems unlikely that these 
clusters could provide the qualitatively different type of 
data required to resolve the questions left outstanding 
in Sec. V. An alternative, but still non-perturbative and 
predominantly unbiased, approach would be the use of 
variational wave functions, either formulated generally 
or in the more speciflc projected wave function technique 
which leads to different types of flux phasci ^"^'^^ Adapt- 
ing this type of treatment to the coupled spin and orbital 
sectors without undue approximation remains a technical 
challenge. 

Within the realm of effective models which could be 
obtained by simpliflcation of the ground-state manifold, 
we cite only the possibility motivated by the current re- 
sults of constructing dimer models based on (ss/ot) and 
(os/st) dimers. Dimer models^ are in general highly 
simplifled, and there is no systematic procedure for their 
derivation from a realistic Hamiltonian, but they are 
thought to capture the essential physics of certain classes 
of dimerized systems. Because QDM Hamiltonians pro- 
vide exact solutions, and in some cases genuine exam- 
ples of exotica long sought in spin systems, including the 
RVB phase and deconflned spinon excitations, they rep- 
resent a valuable intermediate step in understanding how 
such phenomena may emerge in real systems. Here we 
have found (i) a very strong tendency to dimer forma- 
tion, (ii) a large semi-classical degeneracy of basis states 
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formed from these dimers, and (iii) that resonance pro- 
cesses even at the four-site plaquette scale provide a very 
significant energetic contribution. From the final obser- 
vation alone, a minimal QDM, meaning only exchange of 
parallel dimers of all three directions and on all possible 
plaquette units, would already be expected to contain 
the most significant corrections to the VB energy. At 
this point we emphasize that, because of the change of 
SU(2) orbital sector with lattice direction, our 2D models 
are not close to the SU(4) point where four-site plaque- 
tte formation, and hence very probably a crystallization, 
would be expectedJ^ From the results of Sees. IV and 
V, a rather more likely phase of the QDM would be one 
with complete plaquette resonance through all three col- 
ors, and without breaking of translational symmetry. 

Rigorous proof of a liquid phase, such as that repre- 
sented by an RVB state, is more complex, and as noted 
in Sec. I it requires satisfying both energetic and topo- 
logical criteria. Following the prescription in Ref. [lOj, 
three conditions must be obeyed: (i) a propensity for 
dimer formation, (ii) a highly degenerate manifold of ba- 
sis states from which the RVB ground state may be con- 
structed, and (iii) a mapping of the system to a liquid 
phase of a QDM. Criteria (i) and (ii) match closely the 
labels in the previous paragraph, and both dimer forma- 
tion and high degeneracy have been demonstrated ex- 
tensively here. The energetic part of criterion (iii) also 
appears to be obeyed here: static dimers have an en- 
ergy {V), and allowing their location and orientation 
to change gains more (t). The regime V/t < 1 of the 
triangular-lattice QDM is the RVB phase demonstrated 
in Ref. 2^, whose properties include short-range corre- 
lation functions and gapped, deconfined spinous. This 
mapping also contains the criterion of togological degen- 
eracy, and could in principle be partially circumvented by 
a direct demonstration. However, no suitable numerical 



studies are available of non-simply connected systems, 
and so here we can present only plausibility arguments 
based on the high degeneracy and spatial topology of the 
dimer systems analyzed in Sees. IV and V. It is safe to 
conclude that the threefold-degenerate t2g orbital system 
on the triangular lattice is one of best candidates yet for 
a true spin-orbital RVB phase. 

In closing, spin-orbital models have become a frontier 
of intense current interest for both experimental and the- 
oretical studies of novel magnetic and electronic states 
emerging as a consequence of intrinsic frustration. Our 
model has close parallels to, and yet crucial differences 
from, similar studies of manganites (cubic systems of eg 
orbitals), LiNi02 (triangular, eg), YTiOa and CaVOs 
(cubic, t2g), and many other transition-metal oxides, ap- 
pearing in some respects to be the most frustrated yet 
discussed. One of its key properties, arising from the ex- 
treme (geometrical and interaction-driven) frustration, is 
that ordered states become entirely uncompetitive com- 
pared to the resonance energy gained by maximizing 
quantum (spin and orbital) fluctuations. In the orbital 
sector, the restoration of symmetry by orbital fluctua- 
tions makes the model a strong candidate to display an 
orbital liquid phase. Because this liquid is based on ro- 
bust dimer states, the mechanism for its formation is very 
likely to be spin-orbital RVB physics. 
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